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ABSTRACT 

We calculate the observable properties of the most massive high-redshift galaxies in the hierarchical 
formation scenario where stellar spheroid and supermassive black hole growth are fueled by gas- 
rich mergers. Combining high-resolution hydrodynamical simulations of the hierarchical formation 
of a z ~ 6 quasar, stellar population synthesis models, template AGN spectra, prescriptions for 
interstellar and intergalactic absorption, and the response of modern telescopes, the photometric 
evolution of galaxies destined to host z ~ 6 quasars are modeled at redshifts z ~ 4 — 14. These 
massive galaxies, with enormous stellar masses of M* ~ 10 11 — 10 12 M Q and star formation rates of 
SFR ~ 10 3 — 10 4 M Q yr _1 at z > 7, satisfy a variety of photometric selection criteria based on Lyman- 
break techniques including V-band dropouts at z > 5, i-band dropouts at z > 6, and z-band dropouts 
at z > 7. The observability of the most massive high-rcdshift galaxies is assessed and compared with 
a wide range of existing and proposed photometric surveys including SDSS, GOODS/HUDF, NO AO 
WDFS, UKIDSS, the IRAC Shallow Survey, Ultra- VISTA, DUNE, Pan-STARRS, LSST, and SNAP. 
Massive stellar spheroids descended from z ~ 6 quasars will likely be detected at z ~ 4 by existing 
surveys, but owing to their low number densities the discovery of quasar progenitor galaxies at z > 7 
will likely require future surveys of large portions of the sky (> 0.5%) at wavelengths A > 1/zm. The 
detection of rare, star-bursting, massive galaxies at redshifts z > 6 would provide support for the 
hierarchical formation of the earliest quasars and characterize the primitive star-formation histories 
of the most luminous elliptical galaxies. 
Subject headings: galaxies: formation - galaxies: evolution 



1. INTRODUCTION 

Supermassive black holes (SMBHs) serve as the rela- 
tivistic engines of quasars (|Lvnden-Belllll969h and the 
existence of quasars at z ~ 6 suggests the rapid for- 
mation of SMBHs with mas s es Mm ~ 10 9 M Q (e.g., 
iFan et alJ l2000l l2?)0lal l200l l20f)i EQQI). Examining 
the consistency of rapid SMBH growth and the forma- 
tion of cold dark matter (CDM) structures at high- 
rcdshifts has been of long stand ing theoretical inter- 
est (e.g., Efstathiou fc Reesl 119881 ). as a confluence of 
the deep potential well of a rare, massive dark matter 
halo, an abundant supply of baryons, and an efficient 
SMBH fueling mechanism is likely needed to grow enor- 
mous SMBHs in the available « 900 Myrs before red- 
shift z ~ 6. The existence of quasars at z ~ 6 can be 
considered a posteriori evidence that our currently dark 
energy-dominated (ACDM) universe can muster this con- 
fluence, but little else is known observationally about the 
associated formation of quasar host galaxies at such high 
redshifts. The focus of this paper is to assist in the de- 
tection and characterization of massive galaxies destined 
to host z ~ 6 quasars by calculating their photometric 
evolution in the scenario where ACDM structure forma- 
tion fuels SMBHs and grows stellar spheroids through 
gas- rich mergers. 
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The formation of stellar spheroids through mergers 
of disk galaxies has b een a long-established vehicle for 
galax y evolution (e.g.. iToomre fc Toomrd[r972t iToomrd 
Il977| ). If the disk galaxies contain gas before the merger, 
gravitational torques can enable gas to loose angular mo- 
mentum, flow efficiently to the ce nters of the merging sys- 
tems, and result in a starburst (Larson & Tinslcy 197! 



Noguchill9 88: Hernquist 1989; Barne s fc Hernquistlll991l 



19961 : iMihos fc Hernauis t 1994, 1996). The efficient fu- 



eling of starbursts by mergers may explain ultralumi- 
nous infrared galaxies CULIRGS. e.g. ISoi fcr et al. 11984 
Young et al.lll98^ lSoifer"et^ll98V7t fSanders fc Mirabel 
1996h as an intermediate phase where the high in- 
frared luminosities are powered by star formation 
before gas accretion can grow central supermassive 
black holes t o be massive enough t o become lumi- 
nous quasars (|Sanders et al.l Il988bllal: iHernquistJ 119891 : 
ISanders fc Mirabell 119961 : iGenzel etaflll998t ). This evo- 
lutionary picture would paint starbursts, ULIRGS, and 
quasars as stages in a structure formation process that 
builds spheroid and SMBH mass contemporaneously 
through galaxy mergers. 

The observed correlations between SMBH mass 
and stel lar spheroid mass (the Mbh~ Mbuige rela- 
tion, e.g. Kormendv & Richstone 1995: Magorrian et al. 
[19981 iMarconi et all 1200.1 lHaring fc R.ixll2004h . SMBH 
mass and stellar velocity di s persio n (the Mbh~ g_ 
relation. iFerrarese fc MerrittJ 120001: iGebhardt et all 
120001 : iTremaine et all l2002D.and the bl ack hole fun- 
damental plane (Hopkin s et alJ l2007af ) have been 
widely interpreted as evidence for the coeval or cou- 
pled formation of SMBHs and stellar spheroids (e.g. 
Silk fc Reesl Il998t iFabiad Il999t iKauffmann fc Haehneltj 
20001 : iVolonteri et all 120031 : iWvithe fc Loebl 120031: iKmgl 
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12001 iGranato et all [200l iHopkins et all l2006allbT l. 

Hydrodynamical simulations of the formation of 
SM BHs and their host spheroids in gas-rich merg- 
ers (|Di Matteo et al J 120051 : IRobertson et al.ll2006b[ ) that 
includ e a prescription for SMBH growth and feed- 
back (jSpringel et al.l l2005bf l generically predict that 
such mergers undergo starbursts, consistent with ear- 
lier simu lations that did not include SMBH model- 
ing ( e.g., iMihos fc Hernquistl 11994 iBarnes k. HernquisH 
1996), and that the spheroid forms primarily before 
the growing SMBHs reach their final mass. Once 
the rapidly growing SMBHs become massive enough to 
drive energetic o utflows, the SMBH fee d back truncates 
star f ormation (|Di Matteo et al.l l2005t ISpringel et al.l 
l2005aj) and becomes visible as a luminous quasar as 
the out flow removes obscuri ng material at the galaxy 
center (jHopkins et al.1 120051 ). Ultimately, the end 
state of this dissipative process is a spheroidal rem- 
nant that close ly matches observed e lliptical galaxy 
struct ural (e.g.. IRobertson et al.ll2006a[ ) and kinematic 
(e.g., ICox et al.l l2006bT ) properties a nd satisfies the 
observed SMBH-spheroid correlations (|Di Matteo et al.l 
l2005URobertson et al.ll2006bl I Hopkins et al.ll2007ah . The 
general consistency of the evolutionary sequence as de- 
tailed by these simulations with available observational 
data at redshifts z < 3 has been discussed elsewhere (e.g., 
IHonkins et aT1l27)06~ffl c!). 

While a developmental connection between stellar 
spheroids and the SMBHs they host has observational 
support in the nearby universe, the properties of z ~ 6 
quasars may complicate this simp le coeval picture. The 
Sloan Digital Sky Survey (SDSS; lYork et al.lf2000h was 
used in concert wi th Two Micron All Sky Survey J-band 
imaging (2MASS: ISkrutskie et al"]ll997h to photometri- 
cally identify the quasar SDSS J1148+5251, sp ectroscop- 
ically confirmed to reside at redshift z = 6.42 (|Fan et al.l 
120031) . Estimates of the SMBH mass from Mgll and 
CIV line widths for this quasar are A /bh ~2-6x 10 9 
([Willott et al.1 120031 : iBarth et al.ll2003D . suggesting that 
the SMBH is near its Eddington limit. 

Radio observations of SDSS Jl 148+5251 infer the pres- 
ence of molecular hydr ogen from CO(3-2) emission, with 
M m = 2.2 x 1O 1O M (I Walter et al.1 12003). The dynam- 
ical mass for the quasar host Mg vn ~ 4.5 x 1Q 10 M(d, es - 
timated from CO line widths ([Walter et al.ll2003ll200l . 
suggests that the molecular gas and SMBH, rather than 
some dominant stellar component, supply the majority 
of the local potential within ~ 2.5kpc. However, radia- 
tive transfer calculations of CO line emission in quasars 
show that owing to projection effects the small observed 
line widths could be consistent with presence of a mas- 
sive bulge in the quasar h ost galaxy, but such ob served 
orientations are unlikely (|Naravanan et all I2006D. The 
large IR luminosity L m k OO 1 ^ (iBertoldi et apl2003l) 
and radio continuum imaging (|Carilli et al.ll2004D imply 
a star formation rate of ~ 3, 000M©yr _1 if the emission 
is powered by reprocessed starlight, and the presence of 
dust and metals may s uggest a prior epo c h of star for- 
mation in the system fBarth et al.l 12003]: | Walter et al.l 
120031: IBertoldi et al.ll2003l: ICarilli et al.ll2004h . The gen- 
eral uncertainty about the presence of a massive stellar 
component in z ~ 6 quasars then poses two questions: 
if the progenitors and descendants of z ~ 6 quasars do 
contain a substantial stellar component what would their 



properties be and how might one detect them? 

The detection of high-redshift galaxy populations at 
redshifts z > 3 has greatly developed over the last fif- 
teen years. The adoption of photometric techniques 
to identify the Lyman break in galaxy spectra (e.g., 
ISteidel fc Hamilto'nllT992l . [1991 ISteidel et al.lll996D . ow- 
ing to intergalactic medium (IGM) absorption (e.g., 
Madaulll995l). the Balmer and 4000A spectral break (e.g., 
Franx et a l. 2003), and narrow-band searches for Lya 
emission Ce.g. iCowie fcHv3fl998l I5u et al.ll2002D has led 
to the discovery of galaxy populations at progressively 
higher redshifts. Multicolor Hubble Space Telescope 
(HST) observations of the Hubble Deep Field (HDF; 
Willi ams et al.1 119961 ) allowed for the efficient identifi- 
cation of large numbers of £?435-band ( F435W filter) 
dropo ut candidates at 3.5 < z < 4.5 (|Madau et al.l 
19961) and redder, mo re distant dropout candidates (e.g., 



Lanzetta et al.|[l996h including a V^o6-band (F606W fil- 



ter) dropout spectroscop ically verified to reside at z 
5.34 (ISpinrad et al.1 [19981) . 

Subsequently, photometric dropout techniques based 
on the Lyman spectral break were used to identify high- 
rcdshift galax y populations in the Hubble Ultra Deep 
Field (HUDF; iBeckwith et~aLll2006l) wit h the Advanced 
Camera for Surveys (ACS: iFord eTaTl Il998h and the 
Near Infrared Camera and Mu lti- Object Spectromete r 
(NICMOS) Ultra Deep Field ([Thompson et al.l f2005h . 
and in the wide r Great Obseryatories Origins Deep 
Survey (GOODS: iDickinson et afll2003l : IGiavalisco et all 
l2004al) . As with the HDF data, color se lections based on 
B435- a nd V606-dropouts were used (e.g. . IGiavalisco et al.l 
2004b), but large samples based on i 775 -band (F775W 
filter) and even z 850 -band (F850LP filter) dropouts 
have been identified i n GOODS and the H UDF b y 
sever a l groups (e.g., IStanwav et alj [20031 l2004blla[ . 
20051: lYalTetall [2003, 200 41: lYan fcW indhorst T20cT 



Yan et al.l 2006: Bunker et al.ll20 03. 2004; Bouwcn s et al 
20031 12004 iBouwens fc Illingworthl 120061: iBremer et al 



2004; 



Dickinson et al 

— TT — 



2004 IStiavelli et al.1 12004 



Evles et al.l 120061: iDunlop et al.l 20061: Overzier et al. 
1 2006 1 : iLee et al.l l2006rjReddv et al.l 120061: iVanzella et all 
120061: iVerma et all 12001 ! The GOODS optical and 
near-infrared dat a sets have been com bined with Spitzer 
Space Telescope (| Werner et al.l 12004) observations us ing 
the Infrared Array Camera (IRAC: iFazio et al.ir2004h to 
help constrain longer wavelength portions of the spectral 
energy distributio ns (SEDs) of high-re d shift galaxy- 
cand i dates (e.g. , iMobasher et al.l 120051: lEeami et al. 
2005 ] lYaTetatl 120051: lEvles et al.l \20Wt iStark et al. 
2006^ iPapovich et al.1 120061 : lLabbe et al.ll2005L 120061 . see 
also iBarmbv et al.l (|2004l ) ) . 

Parallel efforts have been made with ground-based ob- 
servations, sometimes in concert with HS T or Spitzer 
data sets, including opt i cal break searches (iFranx et al.l 

^ 20041) . z > 



data sets, including opt ic 
12001 iDaddi et al.l l2003: 



van Dokkum et al 



5 Lyman-break galaxy (LBG) sea rches using a va- 
riety of techni ques including R Iz ([Leh nert & B remerl 
200J, BVRiz (lOuchi et al.ll2004h. BzK , (Grazian et all 
20061) . and zSZJHK ([Richard et all |200tJ color ~ 
lections, and very-high redshift (z > 5) na rrow band 
and s pectroscopic Lyq-em itter search es (e.g.. lEllis et al.l 
1 20011: iKodaira et al.ll2003t iKurk et al.ll2004lSantosll2004 
buchi et all 120051 : [chary et all 120051: iTaniguchi et al j 
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2005tlStern et al.ll2005bl ; [lve et al.ll2006HKashikawa et all 



20061 ). Collectively, these observational efforts and oth- 
ers have established the existence of a substantial galaxy 
population at z > 5. A brief review of the known high- 
redshi ft galaxy population can be found in lHu fe Cowid 
(I20M ). 

Theoretical calculations of the photometric properties 
of galaxies have been explored previously by hydrody- 
namical simulations of cosmological s tructure formation 
Katz et all H99I 



at redshifts z < 3 



2003; 



> 



IDave et all [1999; 
2005aUbh. z ~ 3 ([Weinberg et al.) booa 
2004L 120061) . z > , 4 dHarford fc GnedirJ 
Finlator et all 120061: iNight et all 120061 and 



Nagaminc et al. 



Nagamine et al 



. 6 dBarton et all 12004 IDave et all 120061) . and 
in semi- analytic calcu l ations (e.g.. iBaugh eHD H99l 
iKauffmann et all Il999t iSomerville et all 120011 ) These 
calculations show that simulated galaxies, typically with 
maximum star formation rates of SFR < lOOM yr _1 , 
can satisfy Lyman-break color selection techniques. 
These simulated galaxy samples broadly agree with the 
observed properties of high-redshift galaxies, but given 
the simulated comoving volumes (< 0.001/?r 3 Gpc 3 ), the 
galaxies examined to date are orders of magnitudes less 
massive than z ~ 6 quasar host halos that have comov- 
ing number densities of n ~ lGpc~ 3 (jFan et al.ll2003[ ). 
Furthermore, these cosmological simulations have not in- 
cluded prescriptions for energetic feedback from SMBH 
growth that have been demonstrated to influence the 
photometric properties of elliptical remnants formed in 
hydrody namical simulations of isolated disk galaxy merg- 
ers (e.g.. ISpringel et al"1l2005af) . 

To calculate the photometric evolution of z ~ 6 quasar 
host galaxies, we use the simulations of iLi et all (2006) 
that combine merger trees measured from cosmologi- 
cal N-body simulations of a l/i~ 3 Gpc 3 comoving volume 
with hydrodynamical simulations of galaxy mergers that 
include gas cooling, star formation and supernovae, and a 
prescription for feedback from growing SMBHs to track 
the hierarchical formation of a massive galaxy at red- 
shifts z < 14. T hese simulations are described in detail in 
ILi et al.l (|2006l ) and summarized in [J2] The galaxy hosts 
a quasar at z ~ 6.5 that is similar in bolometri c luminos- 
ity to the z = 6.42 quasar SDSS J1148+5251 (|Fan et all 
2003). The photometric evolution of the host galaxy 
is modeled before, during, and after the quasar phase 
by applying population synthesis models, template AGN 
spectra, a prescription for dust attenuation and redden- 
ing, intergalactic medium absorption, and telescope re- 
sponses, and is described in In SJU the photometric 
evolution of the host galaxy is compared with a variety 
of optical and near-infrared photometric selection crite- 
ria from the literature. The observability of high-redshift 
quasar host galaxies is assessed in Sj5] We provide a 
brief discussion in Sj6] and summarize and conclude in [J71 
We adopt a flat ACDM cosmology with matter density 
Ojvl = 0.3, dark energy density Oa = 0.7, baryon density 
Ob = 0.04, spectral index n = 1, and root-mean-squared 
fluctuations on 8/i _1 Mpc scales of a% = 0.9. Unless oth- 
erwis e noted, we report magnitudes on the AB scale (|Okd 
fl97l . 

2. SIMULATIONS 

We use the lLi et al.l ()2006f ) hydrodynamical simulation 
of the formation of a z ~ 6 quasar as input for our cal- 



culation of the photometric evolution of the most mas- 
sive quasar host galaxies. The simulation is described 
in detail by ILi et al.l (|2006l ) , and other sp e cifics o f the 
modeling can b e found in ISpringel et al.l (|2005bf ) and 
iRobertson et al.l j2006b) , but a brief summary of the sim- 
ulation follows. 

The rarity of z ~ 6 quasars requires the simulation 
of very large volumes to capture the high-sigma fluctu- 
ations that give rise to early- forming massive galaxies. 
To accomplish this, a 400 3 particle N-body simulation of 
a l/i~ 3 Gpc 3 volume of a cosmology consistent with the 
first-year results of the Wilkinson Microwave Anisotropy 
Pro be (WMAP1, M = 0.3, A = 0.7, fl h = 0.04, a& = 
0.9. lSDergei"eTall 20031 was performed. This simulation, 
with dark matter particle mass jtidm = 1-3 X 1O 12 /i _1 M 
and force resolution e = 125/i kpc, was evolved from 
z ~ 30 to z ~ 0. Dark matter particles bound to 
the most massive halo in the volume, with Mhaio = 
3.6 x lO^/i^M© at z = 0, are id entified at z » 6 . 
The Gpc 3 volum e is then resimulated (jPower et alJl2003l : 
iGao et all 120051 ) from z = 69 to z = with 350 3 dark 
matter particles, with 340 3 high-resolution dark matter 
particles (M dm = 2.8 x 1O 8 /i _1 M , e = 5/i _1 kpc) dis- 
tributed in the ~ 1.25 x 10~ 4 /i~ 3 Gpc 3 Lagrangian vol- 
ume that comprised the most massive halo at z ~ 6 and 
the remaining particles properly placed throughout the 
rest of the Gpc 3 volume to retain the large scale tidal 
field. 

The ACDM simulation is used to follow each merger 
in the hierarchical formation history of the halo between 
z = 14.4 and z = 6. The merger history is comprised of 7 
major mergers (with mass ratios less than 5 : 1) between 
z = 14.4 and z = 8.5 contributing to the 5. 4x 10 12 /i _1 i\/ Q 
halo mass at z ~ 6. To model the formation of the 
baryonic component of the quasar host halo, each event 
in the merger tree is resimulated using the smoothed 
particle hydrodynami cs (SPH) / N-body code GAD- 
GET2 (lSpringelll2005T) . The merging systems are rep- 
resented with hydrody namical models of disk galaxies 
(|Springel et al.l l2005b| ) with structural p roperties ini- 
tialized according to the iMo et all (|1998f ) dissipational 
galaxy formation forma lism and scaled appropri ately to 
high-redshift following Ro bertson et al.l (|2006bf ). The 
scalings simply fold the redshi ft- and mass-depen dent 
dark matter concentrations (e.g. lBullock et aD2001l ) and 
the redshift-depend ent Hubble parameter H{z) through 
the lMo et all (|1998f) formalism, resulting in smaller, more 
centrally concentrated disks at higher redshifts. The or- 
bi tal parameters of each progenitor galaxy (see Table 1 
of ILi et all ( (20061) ) are set such that the initial separation 
of the galaxies is the current host halo virial radius at the 
time of the merger and the pericentric passage distance 
is half the progenitor disk scale-length i? per i = 0.5i?d- 

The merging disk galaxies range in virial mass from the 
initial progenitor's M vir = 6.3 x lO lo /i _1 M at z = 14.4 
to the last merging galaxy's M vlr = 2.08 x lO 12 /i _1 M . 
Each progenitor disk has a mass mj = 0.15M v i r , with gas 
fractions of / gas = 1 at redshifts z > 10 and / gas = 0.9 
at redshifts 8 < z < 10. For the numerical parame- 
ters of the progenitor galaxy models, a fixed dark matter 
particle mass (Md m = 1.1 x 1O 7 /i _1 M ) and baryonic 
particle mass (M gas = M star = 2.2 x Vfh^Mo) is used. 
The gravitational softening for the dark matter particles 
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(edm = 60/i _1 pc) is set twice as large as for the baryonic 
par ticles (e gas . stars = 30 / i~ 1 pc ) . Star formation is treated 
in a lMcKee fc Ostriken (|1977t ) multiphase picture for the 
interstellar medium (ISM), implem e nted as a numerical 
algorithm by ISpringel fc Hernquistl ([2003) . In this mul- 
tiphase ISM model, dense, star-forming gas particles are 
treated as a collection of cold molecular clouds embed- 
ded in a tenuous phase heated by supernovae. The star- 
formation rate i s calibrated to reproduce the Sc hmidt- 
Kennicutt Law (jSchmidt] [19591 : iKennicuttl 1 19981) . with 
the star forming gas mass calculated from the cold phase 
density. Stellar particles are spawned probabilistically 
from star- forming gas particles. Each stellar population 
is formed with half the mass of its parent gas particle and 
recycles 2% of its mass as metals throug h supernovae (for 
details . seelSprmgel fc Hern auist 2003)). 

The Hi et all (|2006) simulation includes a prescrip- 
tion for SMBH accre tion and feedback (for details, see 
ISpringel et alj l2005bf h Each simulated galaxy is ini- 
tialized with a SMBH "sink" particle that is allowed 
to spherically accrete through a iBondil (|1952D prescrip- 
tion. The accretion rate is determined from the gas 
temperature and density near the SMBH particle, and 
is limited to the Eddington rate. A fraction e r = 0.1 
of the accreted rest mass energy released as radiation 
and a small fraction of this radiated energy (5%) is cou- 
pled to the surrounding gas as thermal feedback. The 
strength of the feedback is calibrated to reproduce the 
loc al Mbh~ Q" relation in i solated mergers of disk galax- 
ies (|Di Matteo et al.ll2005| ). The mass of the SMBH seed 
in each galaxy is determined by assuming a single 200 A/© 
black hole formed in each halo at z = 30, in accord with 
expectations from Population III star formation (e.g. 
Abel et al.l 120021 : iBromm fc Larson! 12001 lYoshida et all 



20031 120061 ). and then experienced Eddington-limited 



growth until its host galaxy merged into the quasar host 
halo. The SMBH seeds initialized in this manner range 
in mass from Af B H = 1.5 x 10 4 /i -1 Af Q at z = 14.4 to 
M BH = 8.92 x 10 6 h- 1 M Q at z = 8.5. As dynamical 
friction operates on the merging galaxies during the sim- 
ulation, the central regions of each system eventually co- 
alesce. Once the SMBH particles approach to within the 
gravitational force resolution they are merged. 

2.1. Mergers, Star Formation Rate, and Active Galactic 
Nuclei Luminosity 

The seven major mergers between eight galaxies oc- 
cur during five redshift eras of the evolution: the first 
merger between two galaxies (M v ; ri i = 5.3 x 10 10 h~ 1 MQ, 
M vir , 2 = 6.3 x I0 w h- 1 M Q ) begins at z = 14.4, followed 
by two additional mergers (M v i r ^ 3 = 1.5 x 1O 11 /i _1 M , 
M vir , 4 = 1.77 x 1O u /i -1 M ) at 13.0 > z > 10.5, single 
mergers at 10.5 > z > 9.4 (M vh , 5 = 4.91 x KF/i^M©) 
and at 9.4 > z > 8.5 (Af viri6 = 7.96 x W 11 h- 1 M Q ), 
and two final mergers (Af viri 7 = 1.6 x 10 12 h~ 1 M (3 , 
M vir , 8 = 2.08 x 10 12 /i _1 M o ) at z < 8.5. The fresh gas 
supply introduced by each merger, combined with the 
high density environment and violent tidal interactions, 
leads to peaks in the star formation rate (SFR, Fig. [TJ 
upper panel) that temporally correlate with the merger 
events. The SFR peaks near the z ~ 8.5 mergers with 
SFR > lO 4 Af yr -1 , but frequently exceeds lO^Qyr" 1 
over the redshift range 8 < z < 13. 
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Fig. 1. — Star formation rate (SFR) and supermassive black 
hole (SMBH) accretion rate as a function of redshift in the 
iLi et al] ((2006) simulation of the hierarchical formation of a 
z ~ 6 quasar. The upper panel shows the SFR determined 
from the density of cold gas over the course of eight galaxy 
mergers in the build-up to quasar activity at z ~ 6.5. The 
lower panel shows the accretion rate Mbh on to SMBHs in 
the system, determined through a lBondil (|1952| ) spherical ac- 
cretion model that depends on t he properties of gas ne ar the 
SMBH particles (for details, see ISpringel et all [2005bl ) . The 
unattenuated bolometric luminosity of the active galactic nu- 
clei (AGN) fueled by SMBH accretion can be calculated from 

panel) also influences the observable properties of the 
system. The unattenuated bolometric luminosity Lboi of 
the active galactic nuclei (AGN) powered by the grow- 
ing SMBHs is determined by the accretion rate onto the 
SMBHs through the equation 



L 



bol 



e r M B HC 



1-15 x ID 1 -' 



AL 



BH 



M Q yr" 



Lr, 



where e r = 0.1 is the radiative efficiency of the accretion 
model we adopt and c is the speed of light. 

3. PHOTOMETRIC MODELING 

Modeling the observed spectral energy distribution 
(SED) of the quasar host galaxy involves calculating 
the SEDs of each individual stellar population and 
AGN present in the system and applying a wavelength- 
dependent dust correction to the calculated SEDs that 
accounts for internal interstellar absorption. The emis- 
sion from the galaxy, attenuated by ISM reddening, must 
then be corrected for both the IGM opacity owing to Ly- 
man series absorption and redshifting of the SED owing 
to cosmological expansion. The measured photometry of 
the object is then modified by the combined wavelength- 
dependent transmission of the filter and telescope assem- 
bly as well as the wavelength-dependent quantum effi- 
ciency of the detector. Below we describe how we model 
each contribution to the observable photometric proper- 
ties of the massive quasar host galaxy. 

3.1. Stellar Population Synthesis 
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The hydrodynamical simulation follows the formation 
of individual stellar populations formed from dense gas, 
as described in fJ2] Each stellar particle represents several 
million solar masses of stars, so an initial mass function 
(IMF) must be selected to represent distribution of stellar 
masses formed when each stellar particle is spawned. The 
IMF also influences the spectral evolution of each stel- 
lar population throug h its time-dependent mass-to- light 
ratio. We adopt the iKroupal (|2001l 12002ft IMF, which 
describes the number N of stars of a given mass m as 
a power-law dN/dm oc mT s . The Kro upa IMF shares 
the same power- law index s w 2.3 as the ISalpeterl (|1955l ) 
IMF for m > 0.5M Q but flattens to s = 1.3 at lower 
masses. The mass-to-light ratio o f the Kr o upa I MF is 
correspondingly lower than for the ISalpet er (1955) IMF, 
by a factor of ~ 2, when normalized over the mass range 
O.IMq < m < 100M Q . 

The metallicity and age of each stellar particle is 
tracked in the simulation, allowing for the use of 
stellar population synthesis models to calculate an 
evolving SEP f o r eac h particle. We use S t arbur st99 
(jLeitherer et all 119991 : I Vazquez fc Leithererl l2005ft to 
model the spectra from each population, given the mass, 
age, metallicity, and assumed Kroupa IMF. The ver- 
sion of Starburst99 we use includes the P adova stellar 
evolution tracks (e.g., iFagotto et alJ Il994ft to improve 
the performance of Starburst99 fo r older (t > O.lGyr) 
stella r populations (for details, see I Vazquez fe Leithererl 
2005). Illustrations of the time-dependent SEDs pro- 
duce d by Starburst99 can be fo und in Figures 4 and 
5 of I Vazquez fc Leithererl (|2005ft . In practice, we pro- 
duce a grid of Starburst99 spectral models at wavelengths 
0.0091/im < A < 160/im for stellar populations over the 
range of ages t = 10 4 — 10 10 yr, spaced logarithmically, 
and for metallicities of Z = 0.02 — 2.5Z@, and then in- 
terpolate the Starburst99 models to produce an SED ap- 
propriate for each stellar particle. 

3.2. Active Galactic Nuclei SEDs 

As discussed in Ej2j the simulations follow accretion 
onto SMBH particles. Given the radiative efficiency we 
adopt (e r = 0.1), we can calculate the bolometric lumi- 
nosity of the AGN associated with each SMBH (Eq.[H). 
Accord ing to the feedback prescription bv lSpringel et alJ 
(|2005bft . a fraction r\ = 0.05 of this bolometric luminos- 
ity is injected as thermal energy in the surrounding gas. 
The remaining e r (l — if) fraction of the radiated energy 
would be visible as nuclear emission in the absence of in- 
tervening material. To model this emission, we combine 
the power-law + x-r ay reflection compon ent SED tem- 
plate for AGN from IMarconi et all (|2004ft with the av- 
erage quasar emission template from lVanden Berk et all 
((2001ft 

The IMarconi et alJ (|2004ft SED template consists of 
a power-law L v oc v a , with a — —0.44 in the range 
0.13/im < A < 1/im set to match the slope inferred 
from stacked SDSS quasar spectra (|Vanden Berk et"aLl 
2001), and a Rayleigh- Jeans fall-off (a = 2) at longer 
wavelengths. The ultraviolet (0.05^m < A < 0.12/xm) 
emission has a powe r -law slope (s = —1.76) deter- 
mined by iTelfer et all (|2002ft based on the HST spectra 
of ~ 200 quasars, while the x-ray emission at energies 
above 1 keV is modeled with a power -law plus PEXRAV 
model (Magdziarz & Zdziarski 1995) reflection compo- 



nent. The normalization of the x-ray component relative 
to the optical pow er-law component is fix ed to reproduce 
the quantity aox (jZamorani et al.lll981ft that relates the 
ultraviolet (A ~ 0.25Atm) to x-r ay (E w 2 keV) flux 
ratio. The IMarconi et all (|2004ft S EP template uses a 
luminosity-dependent aox following lVignali et al.l (2003) 
that leads to a decreasing fraction of the bolometric lumi- 
nosity emitted in x-rays as the AGN luminosity increases. 
The ratio of bolometric luminosity to B-band luminosity 
Lboi/ 'vbLu^b has a corresp onding decrease with bolomet- 
ric luminosity. Figure 3 of lMarconi et al.l (|2004ft presents 
illustrations of their AGN SED template and luminosity- 
d ependent i3-band an d x-ray bolometric corrections. 

IVanden Berk et all (|2001ft produced a composite 
quasar spectrum by stacking the spectra of more than 
2000 SDSS quasars at wavelengths 0.08^im < A < 
8555/im, including systems with r' absolute magnitudes 
of -18.0 > r' > -26.5 over the redshift span 0.044 < z < 
4.789. The continuum of the geometric mean spectrum 
of these quasars can be represented by a broken power- 
law L v oc v s with s = —0.44 over the wavelength range 
0.13/.tm < A < 0.5/im and s = —2.45 at longer wave- 
lengths, and is re presentative of the composite spectrum 
of z ~ 6 quasars (|Fan et alj |2004) . While the short wave- 
length slope is flatter t han the slope s ~ —0. 9 found for 
high-r e dshift quasars ( Schneider et al.l 119911 : iFan et aTl 
2001bt ISchneider et al.ll200l . the difference can be ac- 
counted for by difficulties in deter mining the contin- 
uum. Prominent emission lines in the lVanden Berk et al.l 
(|2001ft template spectrum include Lya, CIV, Mgll, and 
Ha. 

To com bine the IMarconi et all (120041 ) SED template 
with the IVanden Berk et al.1 ( 20011 ) quasar spectrum, 
we replace the opt ical pow er-law component of the 
IMarconi et all (|2004ft with the lVanden Berk et alJ (|2001ft 
spectrum over wavelengths O.ll^m < A < 0.86^m and 
match the Rayleigh-Jeans and UV power-law continua 
to the emission line spectrum. The spectrum is then 
integrated and renormalized to maintain the desired 
bolometric luminosity, and approx imately mainta i ns the 
emission line equivalent widths and IMarconi et all (|2004l ) 
luminosity-dependent bolometric corrections while scal- 
ing the emission line flux with the AGN luminosity. As 
discussed below, the choice to combine these spectral 
templates only influences our results during the evolu- 
tionary phases when the AGN dominates the emission of 
the system, over small redshift ranges n ear z sa 6.5 an d 
z w 6 (see Fig. Q] above and Fig. 11 of iLi et alJl2006ft . 
but allows us to incorporate t he effect of Lyq emission 
on quasar color selection (e.g-. lFan et al.ll2001al ). 

3.3. Interstellar Reddening and Absorption 

The methods for modeling the emission from stellar 
populations ( ij3.1[) and AGN ( £|3.2j) allow for an esti- 
mation of the unattenuated SEDs of simulated galaxies. 
The gas-rich galaxies involved in the mergers that lead 
to quasar activity at z ~ 6 contain substantial columns 
of hydrogen and, owing to the efficient star formation 
in dense gas, large amounts of obscuring dust. The 
wavelength-dependent attenuation from dust reddening 
should therefore be accounted for when modeling the ob- 
served SED of the quasar host galaxy. 

To incorporate dust attenuation in the photometric 
modeling, the hydrogen column density Njj and an esti- 
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mate of the gas metallicity is calculated along the line- 
of- sight to each s t ellar p opulation. We use the method 
of iHopkins et alj (|2005l ) to calculate the gas density 
on a fine grid from the SPH formalism and the par- 
ticle distribution. For dense gas above the threshold 
for star formation a two-p hase decomposition of the 
ISpringel fc Hernquistl (|2003f) multiphase ISM model is 
performed, assuming the hot and cold phases are in pres- 
sure equilibrium, in order to calculate the hot phase col- 
umn density. For each AGN in the system the column 
density and metallicity is measured along 100 lines-of- 
sight to calculate a distribution of optical depths, allow- 
ing for an estimate of the scatter in the internal galactic 
absorption. For discussions on general agreement be- 
tween the bound-free and metal-line absorption calcu- 
lated from this method for estimating Njj columns to 
AG N and models of the cosmic x -ray background, see, 
e.g. . IHopkins et all (|2006all2007bf ) . 

Once the Nu columns and gas metallicities are calcu- 
lated, the visual attenuation Ay owing to dust is calcu- 
lated by scaling the Milk y Way (MW) dust- t o-gas ratio 
A V /N H = 5.35 x 10~ 22 (iBohlin et al J 11978. [Cm3l2OO0h 
by the gas metallicity. A linear scaling with metallic- 
ity, normalized to the MW ratio for solar metallicity gas, 
roughly accounts for the observed dust-to-gas ratios in 
the less metal-rich Large Magellanic Cloud (LMC) and 
Small Magellanic Clo ud (SMC) (e.g., lFitzpatrid3 ll^M 
iRodrigues et al.|[l997l ). The wavelength-dependent dust 
attenuation A\/Ay must also be modeled , with possible 
choic es including the MW dust curve (e.g.. lCardelli et ail 
1989) that includes the strong 0.22/,m absorption fea- 
ture LMC and SMC dust (e.g., |Pelfl99l iGordon etail 
l2003h . and the po ssibly supernoyae-rela ted dust in high- 
rcdshift quasars (iMaiolino et al.ll2004D We adopt the 
dust model advocated by ICalzetti etafl (|1994f ) for mod- 
eling the UV emission fr om starburst galaxies , with the 
updated Ry = 4.05 f rom ICalzetti et all (120001) (see also 
ICalzettilll99"7t ). The ICalzetti et all (|2000D dust law is 
frequently u sed to model LBGs and infer their stella r 
masses (e.g., iPapovich et "all 120011 : IShaplev et a"T1l2005[ L 
an d estimates of unatte nuated star formation rates us- 
ing [ciliittLeFjb (2000) dust corrections have been used 
to successfully pre dict the x-ray and radio properties of 
galaxies at z ~ 2 (jReddv fc Steide]ll200l . Using SMC 
dust would increase the UV absorption of the model spec- 
tra and increase the transition redshift where model spec- 
tra cross boundaries in color-color selections designed to 
be sensitive to the Lyman break, but would not change 
the photometric selection of the model spectra with a 
given color criterion over most of the system's redshift 
evolution. 

3.4. Intergalactic Medium Absorption 

Intervening hydrogen clouds in the intergalactic 
medium (IGM) can greatly attenuate flux via Lyman- 
series abs orption . To a ccount for this IGM opacity, we 
adopt the Madau (1995) model of Poisson-distributed ab- 
sorbers in which Lyman series blanketing causes a "stair- 
case" absorption profile blue- ward of Lya. The spectral 
break introduced by this absorption has been frequently 
used to engineer the photometric selec t ion of high- 
redshi ft galaxies (e.g., iMadau et al.l 119961 : ISteidel et al.l 
Il996t iGiavalisco et al.l [2004bf ) . By z 4.5, the mean 
transmission is only 30% near Lya, but the Poissonian 



nature of the model can lead to substantial variation in 
the optical depth r at intermediate redshifts. At z ~ 3.5 
the ±1(7 fluctuation in optical depth corresponds to 30% 
transmission variation r elative to the unattenuated spec- 
tra (see, e.g., Fig. 3 of lMadaulfl995l ). However, Monte 
Carlo simulations suggest that the bandpass-averaged 
fluctuations in transmis sion are typically sma ller (~ 0.1 
magnitudes, see, e.g., iBershadv et al.lll999D . At high 
redshifts (z > 6) the absorp tion is very strong near Lya 
(e.g., r > 3, see Fig. 2 of iNight et all 12006) and typi- 
cal line-of-sight fluctuations will lead to strong spectral 
breaks. We therefore assume the mean IGM opacity cal- 
culated from the model since we are primarily concerned 
with high-redshifts (z > 5), but we note that the cal- 
culated photometric trajectories could experience either 
small shifts in color-color space position or rate of change 
with redshift owing to line-of-sight opacity fluctuations. 

3.5. Telescope, Filter, and Detector Response 

To calculate the observed photometric properties of 
simulated galaxies, the wavelength-dependent transmis- 
sion of telescope assemblies, filters, and detector quan- 
tum efficiency (QE) must be modeled together to de- 
termine the net transmissivity of a given observational 
facility. 

For the SP S S u' -, g'-, r'-, i'-, and z'-band 
(|Fukugita et al.l 119961 ) transmissivities, we use the 
USNO40 response functions measured by Jim Gymn 6 . 
The calculated U n , G, and R s (jSteidel fc Hamiltodll993h 
magnitudes use the transmissivity measured with the 
KPNO Mosaic CCDs 7 . For comparison with HST ob- 
servations, we calculated the wavelength response curves 
from the combined transmissivity of the mirror and fil- 
ters with the detector QE for F435W, F606W, F775W, 
F814W, and F850LP observations with ACS 8 and 
F110W and F160W observations with NICMOS 9 . For 
the ground-based J-, H-, and i^ s -band transmissivities, 
we us e the calculat e d rela tive response curves for 2MASS 
from iCohen et all (|2003f) . The IRAC throughput re- 
sponse was adopted from data provided by the Spitzer 
Science Center 10 . In practice, over the redshift range 
(z > 4) considered at flux in all filters bluer than F435W 
were always too heavily absorbed to be easily measured 
(> 30 AB magnitudes), so in what follows we focus on 
the observable properties of galaxies in redder filters. 

3.6. Combined SED Modeling 

The final model SEDs calculated for the simulated 
quasar host galaxy combines the stellar population syn- 
thesis models described in tj3.ll the AGN spectral tem- 
plate described in £|3 . 21 ISM attenuation described in 
H3.31 and IGM attenuation described in H'SAl to arrive 
at the observed flux density (adapting the notation of 
IPapovich et 111120011) 



Y^h Lv 0t i(\o,ti, A\ g . i: M+.i, Mbh, 



F u {z,X) = (l + z) 



^D\(z) 



xe 



— tigm(2,A) 



(2) 



6 http:/ /www-star. fnal.gov/ugriz/Filters/response. html 

7 http:/ /www. noao.edu/kpno/mosaic/filters/filters. html 

8 ftp:/ /ftp. stsci.edu/cdbs/cdbs2/comp/acs/ 

9 ftp:/ /ftp. stsci.edu/cdbs/cdbs2/comp/nicmos/ 

10 http:/ /ssc. spitzer. caltech.edu/irac/spectral_response. html 
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Fig. 2. — Schematic diagram of Lya-break dropouts in 
optical and near-infrared filters (F606W, F775W, F850LP, 
F110W, Ff60W, K a , 3.6/xm, 4.5^m, 5.8^m, and 8.0/xm) 
at redshifts z ~ 4 — 14. Cosmological expansion enables 
intergalactic medium (IGM ) absorption owing to hydrogen 
Lyman-series opacity (e.g., Qyiadau 1993) to absorb light at 
progressively longer observed wavelengths. As the Lyman- 
series absorption (~ llOOA, white dashed line) shifts into 
redder filters, galaxies begin to drop out of the correspond- 
ing photometric samples (shaded regions represent a > 20% 
flux decrement in a given band for a flat spectrum). The 
solid lines show the effect of cosmological redshifting on the 
rest-frame wavelength coverage of commonly used filters. 

where the rest-frame emission L VOt i as a function of rest- 
wavelength Ao from each of N sources may depend on 
the its stellar age ti, stellar mass M*^, or SMBH accre- 
tion rate Mbh,i, and the dust attenuation A\ j along 
the line-of-sight to each source. The typical number of 
individual stellar population sources is N > 750, 000 by 
z ~ 6. The total emitted spectrum from the galaxy is at- 
tenuated by the IGM optical depth tigm, and the factors 
of (1 + z) and the luminosity distance Dl(z) account for 
the cosm ology-dependent redsh ift and distance. Again 
following iPapovich et al.l (|200lD . the flux density in each 
filter is calculated as 



f T u F v {z,\)<lv/v 
J T v Avjv 



(3) 



where unit-free T v represents the combined transmissiv- 
ity of the telescope and filter and the detector QE, as 
described in 13.51 The flux density is converted into a 
magnitude on the AB scale (|Okdll97^ ) through the re- 
lation 

m AB = -2.5 log + 23J ' ( 4 ) 

where (F u (z)) is the bandpass-weighted flux density (Eq. 
[3]) and l^Jy = 10~ 29 ergs cm _2 s _1 Hz _1 . 

The machinery described above allows for a thor- 
ough determination of the photometric properties of a 
simulated galaxy and an application of common color 
"dropout" selection techniques based on spectral breaks 
to track the photometric evolution of the quasar host 
galaxy with redshift. For reference, a schematic diagram 
of the relevant redshifts for dropout samples in the ACS 



and NICMOS filters and a characteristic wavelength cov- 
erage for optical and near-IR filters with redshift is pro- 
vided in Figured] 

4. PHOTOMETRIC EVOLUTION 

The photometric evolution of the quasar host galaxy 
is calculated using the methods described in In our 
modeling, we track the luminosity of only the main host 
halo of the quasar during its evolution; the luminosi- 
ties of the merging systems are not followed before they 
enter the host halo virial radius. The merging systems 
are modeled as almost completely gaseous, resulting in 
a smooth photometric evolution as each new system is 
added to the simulation. The variations in the system 
brightness owe to increases in the star formation rate 
during the mergers and the evolution of the combined 
stellar populations of the quasar host galaxy. Below, we 
describe the photometric evolution of the simulated sys- 
tem over the redshifts z — 4 — 14. 

Figure [3] shows the SED and 15 magnitudes calcu- 
lated for the galaxy at six characteristic times [z ~ 12, 
z ~ 10, z ~ 8.5, z ~ 6.5, 2^5, and z ~ 4) dur- 
ing three distinct phases in its evolution. At redshifts 
z > 7, the system rapidly forms stars as the progenitors 
of the quasar host galaxy merge into the system. Each of 
these gas-rich mergers adds enough fuel to induce large 
star formation rates of SFR > 10 3 A/ Q yr _1 . This star- 
burst phase coincides with the formation epoch of the 
stars that later form the quasar host galaxy. Over this 
timespan the stellar mass in the system has grown to 
M* = 1.595 x 1O 12 M , or 96% of the final stellar mass 
at z = 4, with a mass-weighted mean stellar age that in- 
creases from (t) = 124Myr at z = 9 to (t) = 185Myr 
at z = 7.5. The upper panels of Figure [3] show the 
characteristic starburst spectrum of the galaxy during 
this phase, with strong rest-frame UV and blue emis- 
sion. The SED evolves most strongly at long rest-frame 
wavelengths as the mean stellar age of the system in- 
creases, but the short wavelength rest-frame emission is 
remarkably similar during this phase. The luminosity 
of the system increases as the star formation rate and 
stellar mass both increase to z ~ 8. The large optical 
depth owing to Lyman-series absorption causes a near 
step-function break blue-ward of the Lyev line. 

By z ~ 7.5, the quasar host galaxy begins to change 
rapidly. At this time a total gas mass of M gm = 
1.21 x 10 11 M Q remains, which amounts to roughly 7% 
of the total baryonic mass of the system. Of this gas, 
30% has already been shocked, heated by fe edback, or 
is oth erwise too diffuse to form stars (e.g., ICox et aT] 
I2006a|) . The remaining 5% of the star-forming baryons 
is dense enough to support a large star formation rate 
of SFR ~ lOOOAfoyr -1 . However, the starburst epoch 
is clearly over, extinguished primarily by the consump- 
tion of the available fuel rather than by some regula- 
tory feedback mechanism. As the star formation rate 
rapidly declines, the stellar population quickly ages and 
increases to a mean age of (t) = 250Myr by z ~ 7 
and (t) — 407Myr by z ~ 6. While the gas sup- 
ply cannot sustain the formation of a large new pop- 
ulation of stars, it can easily support the exponential 
growth of SMBHs to M B n ~ 10 9 M Q . As the violent 
merging of the system continues to drive the remaining 
dense gas toward the galactic centers a rapid phase of 
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Fig. 3. — Photometric evolution of the quasar host galaxy with redshift. The spectral energy distributions (SEDs), plotted in 
the observed frame, are calculated from stellar and active galactic nuclei template spectra, attenuated for interstellar absorption, 
Lyman-series absorption in the intergalactic medium, and cosmological redshifting and distance (see Ej3j> . At early times (z ~ 
8.5 — 12, upper panels), the galaxy is rapidly star- forming and has ultraviolet-bright rest-frame spectra typical of starburst 
galaxies. By z ~ 6.5 (bottom left panel), the SMBH accretion has peaked, leading to a luminous quasar phase that satisfies 
the Sloan Digital Sky Survey quasar color selection. As the SMBH accretion drops and the massive bulge ages (z < 6, bottom 
middle and right panels), the quasar host galaxy becomes a progressively redder luminous galaxy with a large Balmer/4000A 
break indicative of an aging stellar population. The large bandwidth of some filters can cause an overestimate of the true flux 
density when they straddle the Lya absorption line (e.g., F110W at z ~ 10). 



Eddington-limited SMBH growth begins and eventually 
leads to a luminous quasar phase at z ~ 6.5 after the 
SMBHs particles coalesce. The SED of the galaxy reflects 
the quasar activity, with th e hybrid IVanden Berk et al.l 
(|2001l )- lMarconi et all (|2004fl model AGN spectrum dom- 
inating the blue emission of the galaxy at z « 6.6 (Fig. 
O lower left panel). The massive stellar component of 
the system, with a mean stellar age of ~ 100s Myr, 
contributes substantially to the emission at rest-frame 
wavelengths Ao > 0.3/mi and instills the SED with a 
Balmer/4000A break between the K s - and 3.6/im-bands. 
As feedback from the SMBH heats the gas in the central- 
most regions of the quasar host galaxy, the fuel source 
for SMBH accretion is depleted and the prominence of 
the AGN contribution to the SED declines. By z pa 5.2 
(Fig. [3l lower middle panel) the massive old stellar popu- 
lation with mean stellar age (t) = 588Myr dominates the 
spectrum, with residual SMBH accretion providing less- 
dominant AGN activity. The rest-frame UV and blue 
optical portion of the SED becomes increasingly feeble 
compared with the observed near-IR spectral region pow- 
ered by mature stars. 

Below z ~ 5, the host galaxy continues to age and 
redden during its final, passive evolutionary phase. The 
star formation rate has declined by approximately three 
orders of magnitude from its peak near redshift z ~ 9. 
By z pa 4.3 (Fig. [3j lower right panel) the massive stellar 
component has reached a mean age of (t) = 873Myr, 
and by z pa 4 the mean stellar age will increase to 
(t) = l.OlGyr. The UV and blue optical emission from 
the galaxy has correspondingly decreased, leaving a SED 
with a dominant Balmer/4000A break. The SMBH ac- 
cretion rate has dropped to Mbh ~ 0.2M Q yr _1 , bringing 
the bolometric luminosity of the AGN activity to ~ 10% 



of the bolometric luminosity of the stellar spheroid. If 
the system were to continue to evolve passively to low 
redshift, the galaxy would continue to redden. The gi- 
ant spheroid would eventually become too red to ap- 
pear in LBG samples, such as those defined by U n — G 
(ISteidel & Hamilto n 1993; Stcidcl et al. 1995|) and U 300 - 
B 450 (jMadau et al.l 1996; St eidel et aDl996l ) criteria, but 
would be i dentified in Distan t Red Galaxy (DRG) sam- 
ples (e.g.. iFranx et all 12003( 1 with (J - K s ) = 3.2 and 
(I — K s ) = 5.6 (Vega) at z ~ 3. The quasar descendant 
would eventually obtain aj'-r' = 1.24 color (calculated 
assuming passive evolution to z ~ 0.2) typical of bright- 
est cluster gal axies in the MaxBC G cluster catalogs of 
the SDSS (e.g. lKoester et alJ[200l . 

The SEDs calculated for each of the three phases in 
the photometric evolution of the quasar host galaxy can 
be used to follow the trajectory of the system through 
various color-color spaces. The magnitudes calculated 
according the method described in £|3 - 51 can be com- 
pared with color criteria from the literature to deter- 
mine whether color-selected samples would find mas- 
sive quasar host galaxies. In what follows, we exam- 
ine whether the simulated galaxy satisfies common color 
criteria used in previous observational work. For conve- 
nience, the quasar host galaxy magnitudes in the salient 
observational filters at twenty redshifts z = 14 — 4.25 are 
presented in Table Q] 

4.1. VeQ6~dropout Selection 

The multi-color HST ACS observations in the GOODS 
fields allow for a variety of color-color selection criteria 
to identify LBG samples. The bluest color Lyman-break 
dropout techniques for HST data that overlap with the sim- 
ulation i^d^lrift_cover^e_arethe V606-dropout selections at 
2 ~ 5. Giavalisc o et alj (12004a! ) presented a Veoe-dropout se- 
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TABLE 1 

Photometric Evolution (AB Magnitudes) 



Redshift F435W F606W F775W i' z' F850LP F110W F160W K a [3.6/jm] [4.5/xm] [5.8fjm] [8/xm] 
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Fig. 4. — Calculated Veoe — 4775 and 1775 — zsso colors with 
redshift for a simulated quasar host galaxy spectral energy 
distribution (SED) . Shown is the location of the SED in this 
color-color space when the system SED has significant AGN 
activity (triangles, 5.7 < z < 6.7) or is p assively-evolving 
(diamo nds, z < 5.7) . Also shown ar e the [Giavalisc o et al.l 
|2004al ) GOODS and lYan et aU (|2005l ) HUDF color selection 
criteria for galaxies at z ~ 5. The simulated galaxy satisfies 
these criteria during its passive evolution phase. At higher 
redshifts, during the epoch of AGN activity, the system will 
move into i775-dropout samples (1775 — 2gso ^ 1-3, see Fig. [5J). 



lection for GOODS by the color criterion 



{[(Vm 



1775) > 1-5 
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0.9 x {irr 


V [(V 60 6™ 


«775 


) > 2.0]} 
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{(V&oe- 


«775 


) > 1-2} 




A 




{(«775- 


Z850 


)< 1-3}, 



«850 )] 



(5) 



where V and A represent the logical "OR" and "AND" oper- 
ations, and used galaxies selected in this manner to estimate 
the global star formation rate density at z ~ 5. A similar 



selection criterion, defined as 

{[(Veoe - irrs) > 1.2] A [(1775 - ^sso) < 0.3]} 
V 

{[(V 60 6 - 1775) > 0.4 + 2.67(i 7 75 - z 8 5o)] 

A [0.3 < (i 7 75 - 2850) < 0.6]} 
V 

{[(V606 " «77 5 ) > 2.0] A [(1775 - Z850) > 0.6]} . 



(6) 



was used by lYan et alj |2005) to find z ~ 5 galaxies in the 
HUDF for examination with coincident Spitzer IRAC obser- 
vations. The application of stellar population synthesis mod- 
eling found that galaxies at z ~ 5 in GOODS have typical 
stellar masses of M 4 ~ 1O 1O M , stellar ages of t ~ lGyr, low 
extinction, and a wide range of metallicities. 

Figure [4] shows the calculated Vme — «775 and 1775 — 2850 
colors for the massive host galaxy at redshifts z > 4. The 
system satisfies the Veoe-dropout color criterion over its entire 
passive-evolution phase at redshifts z < 5.7 when the SMBH 
accretion and star formation rates are low and the stellar 
component is passively aging. In Veoe-dropout samples, the 
massive descendants of z ~ 6 quasars could appear as bright 
objects (1775 ~ 23.5) with red K 3 — [3.6^tm] colors (see Fig. [3] 
or Table [IJ. At higher redshifts (z > 5.7), during an epoch of 
AGN activity, the system will move into 1775- and i'-dropout 
samples (1775 — zsso > 1-3, see Fig. [5] and ^4.2|l . 

4.2. inK-dropout Selection 

Earlier in the evolution of the quasar host galaxy, during 
its phase of significant AGN activity at times 5.7 < z < 6.7, 
the IGM absorption has moved into the F775W filter and the 
system becomes red in the 1775 — zsso color. To identify 4775- 
dropou ts at z ~ 6 in the GOODS A CS data, fOiavalisc o et al.l 
( 2004b) and iDickinson et al.l (2004) used a simple criterion 
given by 



(i775 — «85o) > 1-3. 



(7) 



The GOO DS data was al s o use d to define a an i775-dropout 
sample by IStanwav et ail (|2003f ) . using similar criteria 



{^850 < 25.6} 
A 

{(«775 - ^850) > 1-5} . 



(8) 
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Fig. 5. — Calculated 1775 — 2850 color and 2850 magnitude 
with redshift for a simulated quasar host galaxy spectral en- 
ergy distribution (SED). Shown is the location of the SED 
in this color-color space when the system SED is roughly 
starburst-dominated (circles, 2 > 6.7), shows significant AGN 
activity (triangles, 5.7 < 2 < 6.7) or is p assively-evolving 
(diamo nds, 2 < 5.7). Also show n are the fGiavalisc o et al.l 
j2004al ) and lStanwav et all | |2003f ) GOODS color selection cri- 
teria for galaxies at 2 ~ 6. The simulated galaxy satisfies 
these criteria during the epoch of AGN activity and, owing 
to its activity and large stellar mass, is considerably brighter 
than the candidate 2 ~ 6 GOODS galaxies. 

(see also iBunker etP al, 2004). Both samples were used to 
estimate the global star formation rate density at z ~ 6. 

Figure [5] shows the 1775 — zgso color and 2sso magnitude of 
the quasar host galaxy during its evolution. The simulated 
galaxy satisfies these criteria during its phase of AGN activity 
and, owing to its activity and large stellar mass (M* = 1.64 x 
1O 12 M0 at z = 6.5), is considerably brighter (Azsso ~ 3) 
than candidate z ~ 6 GOODS galaxies (zsso ~ 25.5, near the 
IStanwav et al.l (2003) magnitude cut). 

4.3. SDSS Quasar Selection 

The identification of quasars by multi-color selection crite- 
ria has been a common practice for n early twenty years (e.g 
IWarren et al. 1987; Irwin et al. 1991), and has proven espe- 
cially fruitf ul for identifying distant (z ~ 6) q uasars in the 
SDSS (e.g., IFan et al.|[2000l. l2001al. [2001. 12009 ). The revised 
selection criterion used bv lFan et al.l (|2003f ) . given by 

{z < 20.2} 

A 

{(*'-«') >2.2} 
A 

{(2' - Jvcga) < 1.5 + 0.35(i' - z' - 2.2)} , (9) 

incorporates i'-dropouts determined by the SDSS data with 
supplemental J-band photometry to distinguish from objects 
with red z — J colors, such as cool stars. Figure [6] shows the 
z ab — -Tvega and zab — z 'ab colors calculated for the quasar 
host galaxy over its evolution, including the phase of strong 
AGN acti v ity at z ~ 6.5. The simulated galaxy satisfies the 
iFan et all (2003) color criteria during its quasar phase near 
z ~ 6.5 and would be color-selected as a SDSS quasar candi- 
date. However, the system is ~ 2 magn i tudes fainter than the 
z* < 20.2 criterion used by IFan et all (|2003l ) to define their 



Fig. 6. — Calculated z' AB — Jvcga and i' AB — z' AB colors with 
redshift for a simulated quasar host galaxy spectral energy 
distribution (SED). Shown is the location of the SED in this 
color-color space when the system SED is roughly starburst- 
dominated (circles, z > 6.7), has significant AGN activity 
(triangles, 5.7 < z < 6.7), or is passively-evolving (diamonds, 
z < 5.7). Also show n is the SDSS co lor selection criterion for 
z ~ 6 quas ars from | Fan et al.l (|2003l ). The simulated galaxy 
satisfies the IFan et al.l (|2003D color criterion during its quasar 
phase near z ~ 6 and color-selected as a SDSS quasar can- 
didate. However, the system is ~ 2 magnitud es fainter than 
the z < 20.2 criterion used bv lFan et al.l (|2003l ) to define their 
quasar sample. 

quasar sample. The quasar phase in this simulation (z ~ 6.5) 
occurs near the beginning of a broader phase of accretion and 
AGN activity (7 > z > 5.5), before the SMBH mass reaches 
its final mass. Had the peak unobscured SMBH accretion oc- 
curred slightly later (e.g., z ~ 5.8 rather than z ~ 6.5), the 
system would have met the magnitude cut for the IFan et al.l 

(2003) SDSS quasar selection. We note that at z > 6.3 the 
system moves near the L-dwarf locus in the z — J vs. i' — z 
color space, similar to the behavior o f the quasar evolutionary 
track examined bv lFan et al.l |2003) to inform their choice of 
color selection criteria. 

4.4. Zs5o~dropout Selection 

The combined ACS and NICMOS observations of the 
HUDF and GOODS fields have allowed f or a search of z ~ 
7 — 8 g alaxies by looking for 2850-dropouts. lYan fc Windhorst! 

(2004) found three 2850-dropout galaxies in the HUDF but 
concluded that their red Jno — J/160 colors made them 
unlikely to reside at z ~ 7. iBouwens et al.l (|2004f ) and 
iBouwens fc Illingwo rth (2006) used zsso-dropout criteria with 
a blue Jno — -Hi60 color cut to find z ~ 7 galaxies more ro- 
bustly. These color-color criteria include the selection 

{(Z850 - Juo) > 0.8} 
A 

{(zs50 - Jno) > 0.8 + 0.66 (Jno - ffieo)} 
A 

{(Jno - Hieo) < 1-2} (10) 

from Bou wens et al.l (I2004T) . th e "liberal" color selection from 
IBouwens fc Illingworthl (|2006D 

{(2850 " Juo) > 0.8} 
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Fig. 7. — Calculated zsso — Juo and Juo — Hieo col- 
ors with redshift for a simulated quasar host galaxy spec- 
tral energy distribution (SED). Shown is the location of 
the SED in this color-color space when the system SED 
is starburst-dominated (circles, z > 6.7), shows signif- 
icant AGN activity (triangles, 5.7 < z < 6.7), or is 
passively-evolv i ng (di amonds, z < 5. 7). Also shown are the 
IBouwens et al.l lj2004f l flight gray). [Bouwens fc Illingworthl 
(2006)- "liberal" (gray), and IBouwens fc Illingworthl (|2006l V 
"conservative" (dark gray) color selection criteria for galaxies 
at z > 7 in NICMOS observations of the Hubble Ultra Deep 
Field (HUDF). The simulated z ~ 6 quasar host galaxy sat- 
isfies these criteria primarily during its starburst-dominated 
evolution phase at redshifts z ~ 7 — 8, but owing to its ex- 
treme star-formation rates and stellar mass it is typically ~ 6 
F160W-magnitudes brighter than HUDF candidates. 

A 

{(z 85 o - Jno) > 0.8 + 0.4 (Juo - H leo )} 

A 

{(Juo - Hun) < 1-2} , (11) 

and the IBouwens fc I llingworthl (|2006f ) "conservative" selec- 
tion 

{(-Z850 - Juo) > 1-3} 
A 

{ (2850 - Jno) > 1.3 + 0.4 (Juo - Hieo)} 

A 

{(Juo - Hum) < 1-2} . (12) 

These photometric selections find < 5 high-redshift galaxies 
in GOODS, suggesting a rapid decline in the density of lumi- 
nous star-forming galaxies and a decrease of the characteristic 
luminosity of galaxies by « 1 AB magnitude from z > 6 to 
z~7-8. 

Figure [7] shows the zsso — Juo an d Juo — Hieo colors calcu- 
lated for the quasar host galaxy with redshift. The lower- 
redshift epochs that cover the AGN activity and passive- 
evolution phases occur below z ~ 6.7 and largely do not 
satisfy the zsso-dropout color criteria. The starburst phase 
at z > 7, when the luminosity of the system is dominated 
by the young stars generated by the enormous star forma- 
tion rates SFR ~ 10 3 — 10 4 , covers the entire j zsso — Juo 
color extent of the Bouwens fc Illingworthl (|2006l ) conserva- 
tive selection. The starbursting quasar progenitor grows to 
M t w 1.6 x 10 12 Mq in stars with a mass- weighted stellar age 
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Fig. 8. — Calculated Infrared Array Camera (IRAC) colors 
with redshift for a simulated quasar host galaxy spectral en- 
ergy distribution (SED). Shown is the location of the SED 
in this color-color space when the system SED is starburst- 
dominated (circles, z > 6.7), shows significant AGN activity 
(triangles, 5.7 < z < 6.7) or is passively-evolving (diamonds, 
z < 5.7). Also shown are the AG ES active galaxy Deselection 
from lStern et all (|2005bf ) and the lBarmbv et al.l (|2004T ) IRAC 
color selection criteria for LBGs at z ~ 3. The simulated 
galaxy satisfies these criteria during most of i ts evo lution and 
would be detected in both the lBarmbv et al.l (I2004D IRAC ob- 
serva tions and the IRAC Shallow Survey (["Eiscnhar dt et aTl 
I2004f ). but is 3 - 4 AB magnitudes fainter than the AGES 
sub-sample of the IRAC Shallow Survey. 

of (t) — 250Myr at z = 7, powering a remarkable luminos- 
ity ~ 6 F160W-magnitudes brighter than zsso-dropout z ~ 7 
candidates in G OODS. If z ~ 6 quasars form in the manner 
suggested by the lLi et al.l (|2006T ) simulations, a zsso-dropout 
or similar selection criterion would be well-matched to find 
quasar progenitors during the peak of their spheroid forma- 
tion. 

4.5. IRAC-based Selection 

In surveys at redshifts z < 4, the IR properties of both AGN 
and galaxies have been us ed to design color-c olor selection 
criteria using IRAC bands. IStern et al.l (|2005aT ) showed that 
the [3.6pm] — [4.5pm] color of broad-line AGN at 1 < z < 4 
are typically redder than galaxies at z < 2. Their selection, 
given by, 

{([5.8^m] - [8.0/im]) > -0.07} 

A 

{([3.6/itm] - [4.5/H) > -0.157 
+0.188 ([5.8pm] - [8.0/im])} 

A 

{([3.6/H - [4.5pm]) > -2.77 
+2.97 ([5.8/tm] - [8.0/tm])} (13) 

identified 90% of the spectroscopically identified type 1 AGN 
in the AGN and Galax y Evolution Survey (AGES). 

iBarmbv et al.l 1)2004! ) used template SEDs to estimate the 
location of galaxies at z > 2 in the [3.6pm] — [5.8pm] vs. 
[3.6/im] — [4.5pm] and [3.6pm] — [8.0pm] vs. [3.6pm] — [4.5pm] 
color-color spaces. Comparing with a sample of confirmed 
z ~ 3 LBG galaxies, they found that the majority of the 
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LBGs resided in the IRAC color-color regions 
{([3.6/tm] - [4.5/xm]) > -0.15} 



and 



{1.1 > ([3.6/im] - [5.8/im]) > -0-4} 



{([3.6/im] - [4.5/im]) > -0.15} 



{1.1 > ([3.6/im] - [8.0/im]) > -0.5} . 



(14) 



(15) 



iBarmbv et al. (2004) combined these criteria with selections 
based on i? s -band magnitude. At the redshifts we consider, 
the quasar host galaxy always has significant IGM attenua- 
tion of the _R s -band and is significantly redder in R a — [ 3.6/im ] 
color than the galaxies examined by [Barmby e t alj (|2004l ). 
We therefore compare only with their IRAC color-color selec- 
tions. 

Figure [8] shows the calculated IRAC colors for the simu- 
lated q uasar host galaxy at redsh ifts z > 4. The gala xy satis- 
fies the lBarmbv et alj (20041 ) and lStern et alj (|2005bf ) criteria 
during most of its evolution. The [3.6/im] — [4.5/im] color be- 
gins blue at redshifts z > 12 when the star-forming SED has 
its rest-frame near-UV and optical emission in the shortest- 
wavelength IRAC bands. As the redshift decreases and the 
stellar populations age, the [3.6/im] — [4.5/im] color increases 
until the Balmer/4000A break passes through these bands 
and the galaxy becomes very red in [3.6/im] — [4.5/im]. At 
lower redshifts, the system becomes relatively stable in its 
[3.6/im] — [4.5/tm] color for the remainder of the simulation. 
The system has a consistently flat [5.8/ im] — [8.0/im] color , 
and would allow for its selection by the IStern et alj (|2005bD 
criteria. Ho wever, the AGES sub- sample of the IRAC Shal- 
low Survey ( Eiscnhardt et al. 2004) has a limiting magnitude 
~ 3 magnitudes fainter than the complete Shallow Survey in 
the [3.6/im]-band, and i s too shallow t o dete ct the simulated 
quasar host galaxy. The Barmbv et al. (2004) color selections 
involve the [3.6/im]-band and have considerably more evo- 
lution than redder bands since the [3.6/im]-band transitions 
through the Balmer/40 00A break (see j]4.6l below). The large 
color-color area of the IBarmbv et al.l (120041 ) criterion would 
allow for the quasar host galaxy to be selected over almost 
the entire duration of the simulation (redshifts z < 12). All 
three IRAC criteria examined would be convenient methods 
to define color selected samples that could include quasar host 
galaxies at z > 4. 

4.6. High-Redshift 4000A/Balmer Breaks 

The presence of 4000A or Balmer breaks in the observed 
SEDs of galaxies has been cited as ev idence for evolved stellar 
populations in galaxies at z > 2 (e.g. iFranx et alj|2003f) . The 
calculated quasar progenitor star formation history (Fig. [TJ 
suggests that a significant population of stars (M* ~ 1O 11 M0) 
can form in such systems by z ~ 12. As these stars age they 
can contribute a small 4000A / Balmer spectral break in the 
galaxy SED, as is apparent in the z ~ 10 SED in Figure [3] 
and may be detected photometrically in the near-IR. 

Figure [9] shows the calculated IRAC [3.6/im] — [4.5/tm] 
vs. [4.5/im] — [5.8/im] colors for the simulated quasar host 
galaxy at redshifts z > 4. The observable near-IR SED is 
very flat over the majority of the galaxy's evolution, remain- 
ing at ([4.5/im] - [5.8/im]) < 0.2 and ([4.5/im] - [5.8/im]) < 
0.27 at redshifts z < 9. At higher redshifts z ~ 9-11 
when the SED would move into Jno-dropout samples, the 
Balmer/4000A break moves into the [3.6/im]-band and causes 
the [3.6/im] — [4.5/im] color to redder by ~ 0.2 — 0.3 mag- 
nitudes. As the Balmer/4000A break moves through the 




0.2 

[4.5]-[5.B] (AB) 

Fig. 9. — Calculated Infrared Array Camera (IRAC) 
[3.6/im] — [4.5/im] vs. [4.5/im] — [5.8/im] colors with redshift 
for a simulated quasar host galaxy spectral energy distribu- 
tion (SED). Shown is the location of the SED in this color- 
color space when the system SED is starburst-dominated 
(circles, z > 6.7), shows significant AGN activity (trian- 
gles, 5.7 < z < 6.7) or is passively-evolving (diamonds, 
z < 5.7). At z ~ 10, the quasar progenitor has already 
been forming stars for > 200Myr and will display a rest-frame 
Balmer/4000A spectral break. The break could be observed 
as a red [3.6/im] — [4.5/im] color in Jno-dropout samples, anal- 
ogous to the K s — [3.6/im] spectral br eak measured in z g^p- 
band dropouts at lower redshifts (e.g.. iLabbe et aLll2006l ). 

[4.5/im]-band, the evolutionary track of the SED reverses di- 
rection at higher redshifts and the system becomes substan- 
tially bluer in [3.6/im] — [4.5/im] as the filters begin to probe 
the region of the SED sensitive to the large star formation rate 
(at z ~ 12, see the upper left panel of Fig. [3]). In contrast, 
the trajectory of local galaxy SEDs with mature stellar pop- 
ulations in color-spaces sensitive to optical breaks continue 
to redden as the SEDs are shifted to higher re dshifts (e.g., 
J 3 — K color at redshifts z > 2.5, see Figure 1 of IFranx et al.l 
12003ft . The extremum of the [3.6/im] - [4.5/im] at z ~ 9 - 11 
owing to the rest-frame optical break makes the color a use- 
ful indicator for the presence of mature stars in Jno-dropout 
samples. 

A completely analogous behavior can occur in lower red- 
shift samples in bluer bands. SED fits to zgso-dropout galax- 
ies in the HUDF at z ~ 7 show evidence for a rest-frame op - 
tical break in their K — [3.6/im] colors dLabbe et al.1 2006). 
as do z « 6 galaxies l|Yan et al.l 120051 ; lEvles et al.l 120051 ) , 
and p ossibly indicate an ev olved stellar component. Re- 
cently, Wiklind et al. (2007a b) have proposed Balmer break 
selection criteria for post-starburst systems at z ~ 5 using 
VLT/ISAAC near-IR J, H, and K s and IRAC [3.6/im] band 
observations of the GOODS fields. Their selection is defined 
by a union of objects in two regions: a "Region A" in the 
J — K 3 vs. H — [3.6/im] color space, given by 

{(J - Ks) < -1.94 + 3.14 (Ks - [3.6/im])} 

A 

{(J-Ks) > -1.90 + 1.27 (Ks - [3.6/im])} 

A 

{(J - Ks) > 1.71 - 0.82 (K s - [3.6/im])} , (16) 
and a "Region B" in the H — [3.6/im] vs. H — [3.6/im] color 
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Fig. 10. — Calculated near-infrared J, H, K s , and In- 
frared Array Camera (IRAC) [3.6/mi] colors with redshift 
for a simulated quasar host galaxy spectral energy distribu- 
tion (SED). Shown is the location of the SED in this color- 
color space when the system SED is starburst-dominated (cir- 
cles, z > 6.7), shows significant AGN activity (triangles, 
5-7 < z < 6.7) or is passively-evolving (diamonds, z < 5.7). 
Also shown ar e the "Region A " and " Region B" color selection 
criteria from iWiklind et al.l l|2007al lbh designed to identify 
post-starburst systems at z ~ 5 via t he Balmer/4000A brea k. 
The quasar descendant satisfies the IWiklind et a"D l|2007al fbT) 
criteria over the redshift range 4.8 < z < 5.5, during the 
post-starbur st, passively-evolving phase of the system's evo- 
lution. The IWiklind et al.l (|2007al l brl criteria are also briefly 
satisfied by the galaxy at z ~ 6.5 during the height of the 
quasar activity. 



space, given by 



{{H- [3.6/im]) > 1.75} 



{(Ks - [3.6/tm]) > 1.20} . 



(17) 



Figure [10] shows the calculated near-infrared J, H, K 3 , and 
IRAC [3.6/im] colors with redshift for the quas ar host galaxy 
SED in relation to Regions A and B from Wikli nd et al 



2007a 




2007a 





The quasar descendant satisfies the IWiklind et al 
criteria over the redshift range 4.8 < z < 5.5, 
during the post-starburst, passively-evolving phase of the 
system's evolution. The comparison suggests that the 
IWiklind et al] l|2007al lbh criteria would select the descendants 
of z ~ 6 quasars durin g the passively- e volving p hase near red- 
shift z ~ 5. While the lwiklind et all (|2007al lbh selection was 
designed to find massive (M* ~ W 11 Mq) and evolved galaxies 
near z ~ 5, the agreement may be somewhat surprising given 
that the evolution of the simulated quasar host SED in the 
J-, H-, K s -, and [3.6/mi]-bands is qualitatively different than 
the evo lving SED models used to design the IWiklind et al.l 
(|2007al fbT) selection. However, we note that over the range of 
redshifts where the quasar host galaxy would be detectable in 
the GOODS and ISAAC data the si mulated system resides in 
the locus of observed galaxies in the Wiklind et al.l (|200Tal lbh 
sample in the (J — K S )-(H — \ 3.6nm])-(K a — [ 3.6/im]) color 
space. We also note that the IWiklind et al . (2007a.Q) cri- 
teria are briefly satisfied by the galaxy at z ~ 6.5 during 
the height of the quasar activity, but such bright sources can 
be discriminated from inactive high-redshift galaxies through 
spectroscopic observations. 



5. OBSERVABILITY 

The comparisons between the redshift-dependent photo- 
metric properties of massive quasar host galaxies and color 
selection criteria in M. 114.61 demonstrate that the SEDs of 
the most massive high-redshift galaxies likely satisfy a vari- 
ety of existing photometric selection techniques. However, 
these high-redshift galaxies are extremely massive and there- 
fore extraordinarily rare. To determine whether existing or 
future galaxy surveys would include such rare systems in their 
samples, the observability of these massive galaxies must be 
estimated. A straightforward quantification of the observabil- 
ity is to calculate the fractional sky coverage and flux limit 
required to include and detect some defined galaxy sample. 
Below, we use the evolution of the quasar host galaxy SED , 
the ACDM cosmological model, and lPress fc Schechtcr (Il974f ) 
theory to estimate these quantities. 

The comoving volume element dV of a spherical redshift 
shell of thickness dz at comoving distance r(z) and redshift z 
can be written 

' ' (18) 



dV = 4nr (z) — dz 
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c is the speed of light, the redshift-dependent Hubble param- 
eter for a flat universe is 



H(z)=H [n M (l + z) :i + QA 



1/2 



(20) 



and Hq is the Hubble constant today. The number of galax- 
ies diV in this comoving volume element is then simply 
diV = n(z)dV, where n(z) is the redshift-dependent comov- 
ing number density of galaxies that comprise the sample. 
Fo r the form of n(z), we use the mass function proposed 
bvlSheth fc Tormenl |j 1999ft to extend the iPress fc SchechteH 
l| 1974ft theory of halo formation to model ellipsoidal collapse. 
Assuming an overdensity threshold 8 C = 1.686 for halo col- 
lapse, we can define the quantity v — y/a5 c /[cr(M)D(z)] that 
relates the root-mean-squared fluctuations a(M) on a mass 
scale M, the linear growth factor 



D(z) = DqH(z) 



(l + z')dz' 
H 3 (z') 



(21) 



with Do normalized such that D(z = 0) = 1, and 8 C to de- 
scribe the rarity of the density peak that collapsed to form the 
halo of mass M. Here the constant a — 0.70 7 was determined 
from cosmological N-body simulations by ISheth fe Tormenl 
(1999). In terms of the peak rarity i>, the comoving number 
density of galaxies above some given mass Mi ow is then 
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where the parameter A — 0.322 is constrained such that 
J °°(dn/dlnm)dlnm = 1 and p = 0.3 is determined from 
fitting to halo mass functions measured in cosmological sim- 
ulations. If we choose Mi ow to correspond to the redshift- 
dependent virial mass of a quasar host galaxy and assume 
that more massive halos will host brighter galaxies, then Mi ow 
defines a luminosity-selected sample of very massive galaxies 
at redshift z with comoving number density n(z). 

To calculate an actual number N of galaxies, a redshift 
interval zi ow < z < Zhigh must be prescribed to define the 
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comoving volume that hosts the sample as 

N(ziow, zhigh, Miow) = / n{M\ aw ,z)AV 

= 4ttc °( M i..^)yV)^' _ (23) 

For a given selection criterion, such as a Lyman-break 
dropout, zi ow will correspond to the redshift at which galaxy 
SEDs satisfy the color criteria. At zi ow the least massive 
galaxy in the sample will have a magnitude tuab that corre- 
sponds to the brightest possible magnitude limit of a survey 
that could detect galaxies with mass M\ ovr . Typically, for 
dropout criteria this magnitude limit would apply to the red- 
der of the two bands (e.g., i775-band for V606-dropouts). How- 
ever, such a survey would have a very small comoving volume 
(since zi ow = Zhigh in this case) and it is desirable to improve 
the magnitude limit to increase the comoving volume of the 
survey. As the magnitude limit of the survey is improved 
to a new m^Bi Zhigh increases to the epoch when galaxies 
with mass Mi ow first became brighter than m' AB . For a given 
photometric selection which sets zi ow , and survey magnitude 
limit ttiab, the fractional sky coverage needed to detect Nd 
galaxies with mass greater than Mi ow at zi ow is 

/ = N d /N(z lovJ , Z high [OTAB], Mow) (24) 

While this is conceptually straightforward, there are two im- 
mediate complications. First, the correspondence between 
the upper redshift limit Zhigh and the magnitude limit wjab 
depends on the redshift-dependent photometric properties of 
galaxies with mass Mi ow . Second, galaxies with mass Mi ow 
at epoch zi ow are evolving and will have some lower mass at 
higher redshifts. Fortunately, the calculations of the redshift- 
dependent galaxy SED described in |3] provide the desired 
Zhigh — ttiab corresponde nce. For the ma ss accretion history 
the merger-tree from the iLi et al] (|2006f ) simulation is box- 
car averaged with window width of dz = 0.3. The resulting 
smooth mass accretion history increases rapidly from z ~ 14 
to z ~ 9 owing to frequent mergers, and has a shape similar 
to the lWechsler et all (|2002P l form for dark matter halo mass 
accretion histories. We note the observability calculation is 
insensitive to the exact shape of the mass accretion history as 
long as the rapid decline in virial mass at z > 9 is reproduced. 

The relation between fractional sky coverage and survey 
limiting magnitude is plotted in Figure [TT] for a variety of 
color selection criteria explored in [3] including Veoe-dropout 
(Eq. [5] ziow ~ 4.5), i 77 5-dropout (Eq. [7] zi ow ~ 5.7), z 8 so- 
dropout (Eq. [12] zi ow ~ 6.8), and [3.6pm] - [4.5pm] > 0.3 
break (zi ow ~ 9.1) samples. For these c olor selection c riteria, 
the quasar host galaxy virial mass in the lLi et al] (|2006h simu- 
lation increases from Mi ow = 5.22 x 10 12 7i _1 M© at zi ow = 9.1 
to Mi ow = 8.01 x 10 12 /i~ 1 M Q at zi ow ~ 4.5. Beginning with 
the brightest possible magnitude limit, corresponding to the 
luminosity of Afi ow -mass galaxies at zi ow , the necessary frac- 
tional sky coverage declines as the sample comoving volume 
increases with limiting magnitude until Zhigh is reached. At 
larger (fainter) magnitudes, the fractional sky coverage be- 
comes flat. The width of the transition region before the frac- 
tional sky coverage flattens is set by the redshift range over 
which the luminosity declines. For the Vqoq-, 2775-, and zsso- 
dropout samples, IGM absorption sets Zhi g h to approximately 
the redshift when Lya moves into the 1775-, zsso-, and Jno- 
bands, respectively. The Lyman break is dramatic at such 
high redshifts and, when combined with the decrease in dr/dz 
and n(z) with redshift, causes the fractional sky coverage to 
flatten quickly beyond the brightest limiting magnitude. The 
larger transition width for the zsso-dropout selection relative 
to the bluer color selections is set by the comparably large 
bandwidth of the F110W filter. For the [3.6pm] - [4.5pm]- 
break selection, the fractional sky coverage artificially flat- 
tens at the [4.5pm]-band magnitude of the system at z ~ 14 



([4.5pm] ~ 26), where the simulation begins. Since these 
galaxies are bright ([4.5pm] < 24 to z ~ 12), only a relatively 
small fractional sky coverage is needed to detect them. How- 
ever, given the decline of dr/dz and n(z) with redshift, such 
a sample would be dominated by systems at z ~ 4. 

Given the rapid decline of the comoving number density 
n with mass Mi ow at fixed redshift, a useful comparison is 
the fractional sky coverag e needed to ob serve systems that 
are less massive than the ILi et al] l|2006f ) simulated galaxy 
by a factor of ~ 2. Figure [TTJ also shows the require frac- 
tional sky coverage for these systems, assuming their pho- 
tometric evolution is similar to massive galaxies and their 
luminosity scales with their halo mass. Such systems are 
clearly ~ lOx more abundant than galaxies twice their mass. 
These systems would have stellar masses of 7x W 11 M @ , and 
SMBH masses A^bh ~ 10 9 AfQ, and could host quasars at 
z < 6. The virial mass of such galaxies is similar to that 
inferred from the c lustering of quasars a t in the 2dF, SPSS, 
and other surveys dPorciani et al.ll2004l: iPorciani fc Norberd 
20061: IWake et all 120041 ; ICroom et al] 120051 : ICoil et all 12009 : 
Shen et aU l2007|) and measurements of the quasar proxim- 
ity effect ilFaucher-Giguere et al] 120071 : iKim fc Croft] 120071 : 
iGuimaraes et alj|2007[). and arequ alitatively similar to the 
system simulated bv lLi et al] l[2006l) and the large se t of z ~ 6 
merger simulations performed by Robertso n et al] (|2006bf) . 
The detection of the host galaxies of such quasars in Lyman- 
break dropout or Balmer-break samples would also provide 
interesting insights into the formation of massive high-redshift 
galaxies (for a possib le sample of such systems at z ~ 5, see 
iWiklind et aLll2007al lbl). 

For comparison, we also show the necessary sky coverage to 
find objects with an abundance fixed to the comoving num- 
ber density (n ~ 4 x 10~ 9 /i _3 Mpc -3 ) of halos with virial 
mass M v ir > 8 x 10 12 h~ x Mq at z ~ 6.5 if their luminosity 
evolution is similar to the simulated quasar host galaxy (Fig. 
11, dashed-dotted lines). Correspondingly, such systems have 
an abundance similar to massive quasar host galaxies when 
they fall in 2775-dropout samples. For higher-redshift selec- 
tions, this estimate of the massive galaxy abundance is more 
favorable for finding the progenitors of z ~ 6 quasars and is 
comparable to our Mi ow = 0.5M v i r abundance estimates. 

The estimated survey parameters for finding high-redshift 
massive galaxies can be compared with the properties of ex- 
isting and future photometric surveys to gauge the likelihood 
for the most massive high-redshift galaxies to be found in ac- 
tual observational samples. The largest existing survey is the 
SPSS, which has magnitude limits of i' ~ 22.3 and z' ~ 20.8 
<|York et al]|2000h over 8452 deg 2 11 . While the SPSS has suf- 
ficient area to capture the most massive galaxies at redshifts 
as high as z > 12, the bright magnitude limits and lack of 
infrared coverage makes the detection of extremely massive 
high-redshift (z > 6) galaxies in the SPSS implausible. Of 
course, as has been beautifully demonstrated, the rarer sub- 
population of massive galaxies in a bright quasar phase has 
been detected in ~ 1/2 the total SPSS area (e.g., iFan et all 
|2003j) as the remarkable luminosity of the z ~ 6 quasar sam- 
ple and the large sky coverage of the SPSS allows for the 
detection of tens of systems. 

Compared with SPSS, deeper but smaller surveys have 
more advantageous parameters for finding Lyman-break 
dropouts at redshifts z > 4.5. For instance, the Bootes 
field of the National Optical Astronomy Observatory Wide 
Peep Field Survey (Jannuzi fc Pevlll999f ) has an 9.3 deg 2 
field with magnitude limits of / = 26 and J = 21. While 
the J-band coverage is likely too shallow to detect high- 
redshift quasar host galaxies, this area-magnitude combina- 
tion should be sufficient to find the most massive galaxies 

http://www.sdss.org/status/ 
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FlG. 11, Estimated survey parameters required to find z J> 7 quasar progenitors (upper panel) and quasar descendants at redshifts 7 ^ z ^ 4 (middle panel), or 

detect quasar host galaxies in the near infrared (bottom panel). Shown is the fractional sky coverage and minimum AB magnitude limit needed to build VgQg-dropout 
(purple), i 7 7 5 -dropout (blue), zggQ-dropout (green), and [3.6um] — [4.5/im] break (red) samples that include a galaxy more massive than the virial mass M vir (solid line) 
or 0.5M vir (dashed line) of the simulated z ~ 6 quasar host from rTJT~et al.| |2UUbi ) - For comparison, the coverage needed to find objects with an abundance fixed at the 
number density of halos more massive than the simulated galaxy at i ™ 6.5 (dashed-dotted line) is also shown. As the dropout selection moves to redder bands and higher 
redshifts, the comoving volume and redshift interval over which massive galaxies satisfy selection criteria decreases. The comoving number density of massive galaxies, 
calculated using the Bheth & Tormen 199!!) mass function, also declines rapidly at high- redshifts . The combination these effects requires large fractional sky coverage to 
find starbursting quasar progenitors at z > 7 (e.g. zgc;n-dropout or [3. 6um] — [4.5um] break selections). The circles show the parameters of the existing Hubble UDF (1775- 
and z g5 Q-band), GOODS (1775-, 285(T' and [4 . 5rim]-bands) , SDSS (i' - and z'-bands), NOAO Wide Deep Field Survey (I- and ./-bands), UKIRT Deep Sky Survey (Z- and 
J- bands for the Galactic Clusters Survey, J- and iv-bands for the Deep Extra Galactic and Ultradeep Surveys), and IRAC Shallow Survey ( [4. 5um]-band) observations. 
Future observations from the Dark Universe Explorer (R+ I + Z band) and the Ultra-deep Visible and Infrared Telescope for Astronomy survey ( Y~-, ./-, H-, and K-bands) could detect 
z ~ 5 — 6. Even wider future surveys with redder sensitivity, such as the Pan-BTARliS Medium Deep Survey (z- and Y-band), Large Synoptic Survey Telescope observations 
(r' Single Visit), and the SNAP Panoramic and Lcnsing su rveys (error bars) co uld find quasar progenitors at z > 8 if their two reddest-bands reach > 22 AB magnitude 
sensitivity. A James Webb Space Telescope Deep- Wide Survey (Gardner et al. 2006) would need to be considerably wider to find the z > 7 starbursting progenitors of z ~ 6 
quasars. Also shown is the fractional sky coverage to detect z ~ 6 quasar host galaxies at 2 ^ 4 in the [4. 5u.m]-band (orange lines), which only requires a small-area survey 
with a [4. Slim] ~ 20.5 magnitude limit. 



at z ~ 4.5 — 5 in U-dro pout samples. The a ccompany- 
ing IRAC Shallow Survey (jEisenhardt et al.ll2004T ). with 8.5 
deg 2 coverage down to [4.5pim] ~ 20.9 would detect these 
systems as well as higher-redshift /-dropouts in the NOAO 
NWDFS (for galaxie s with photometric redshifts z < 2, see 
iBrodwin et al.l 12006). The United Kingdom Infrared Tele- 
scope Deep Sky Su rvey (UKIDSS, IWarren fc Hewitt! l200l ; 
IWarren et al.|[2007T ) includes a Galaxy Cluster Survey (CCS, 
70.9 deg 2 area to Z ~ 21 and J ~ 20.4), a Deep Extragalactic 
Survey (DXS, 35 deg 2 area to J w 23.2 and K = 22.7), and 
an Ultradeep Survey (UDS, 0.77 deg 2 area to J w 25.7 and 



K w 24.7, see lLawrence et al.l2006l ,. While the UKIDSS CCS 
is too shallow to find the very massive galaxies we model, the 
UKIDSS DXS and UDS could detect such systems in the J- 
and is"-bands at z > 4. The UKIDSS DXS may even have 
enough area to detect massive quasar progenitors at z > 7 
(see also IWarren fc H cwctt 2002), which would be an exciting 
prospect, but the lack of deep coverage at shorter wavelengths 
in the DXS may limit its ability to identify such systems as 
Lyman-break dropouts. A proposed UDS with the Visible 
and Infrared Survey Telescope for Astronomy (VISTA) could 
reach 5-ct AB limiting magnitudes of Y = 26.7, J = 26.6, 
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H = 26.1, and K a = 25.6 over a ~ ldeg 2 field (M. Franx, 
private communication). While the proposed VISTA obser- 
vations would be significantly deeper than the UKDISS UDS 
or the NDWFS, the size of the field could limit the survey to 
finding z ~ 6 quasar host galaxies at z < 5. 

The GOODS observations that cover 365 arcmin 2 reach 

depths of t 7 75 ~ 27.4 and 2 8 so « 26.9 with ACS 

i Giavalisco etH] l2004bl) and [4.5pm] = 24.5 with IRAC 
I Dickinson et all 120031 ) The GOODS V 60 6-dropout sample 
has the necessary sky coverage to include a M* ~ 5 x 1O 11 M0 
stellar mass galaxy at z ~ 4 — 5 with 1775 ~ 25 and detections 
in all redder bands, though we are not aware of any such 
systems in the current GOODS data. Less massive (M* ~ 
1O 11 M0) systems are much more abundant and our estimate 
would predict a handful of such 2 ~ 5 systems in GOODS, 
which may have already been detected (e.g., IStark et al.l 
2006). The HUDF data, with an 11 arcmin 2 area and magni- 
tude limits of 2775 = 29.21 and zsso = 28.67 (Bcckwith et al. 
2006), is simply too narrow to include objects as rare as z ~ 6 
quasar host galaxies. We note briefly that these simulations 
suggest that to detect M* ~ 5 x 10 11 Mp) galax i es in Jno- 
dro pout samples at z > 7 fe.g . iMobasher et al.l 120051 . see 
also lChen fc Marzkd (|2004h and lDunlop et al.l (|2006r n7a sur- 
vey 3-4 orders of magnitude larger than the HUDF would be 
necessary. 

Future surveys with wide area coverage and deep imag- 
ing have a substantially better opportunity to detect the 
starbursting progenitors of z ~ 6 quasars. The Pan- 
STARRS project is planning a 1200 deg 2 Medium Deep Sur- 
vey with magnitude limits of zo.89/jm = 24.7 and Yi.o2/jm = 
23.9 12 , which could find massive quasar progenitors in z- 
band dropout samples at z > 7. The Large Synoptic Survey 
Telescope would make short work of finding the most mas- 
sive high-redshift galaxies if its Y-band Single Visit Depth 
is comparable to its target r = 24.5 depth 13 . The Super- 
novae Acceleration Probe (SNAP) is planning a weak lens- 
ing survey of 300 — 1000 deg 2 to itiab = 28 and a wide-area 
(7000 - 10000 deg 2 ) survey with a depth of m AB = 27 14 . If 
a near-IR HgCdTe device flies with the mission and reaches 
similar sensitivities down to A ~ 1.5/im over the same area, 
SNAP could detect the progenitors of z ~ 6 quasars at red- 
shifts approaching 2 ~ 10. The Dark Universe Explorer mis- 
sion will survey 20, 000 deg 2 with a composite R + 1 + Z-band 
in the red o ptical and will reach a n effective I ~ 24.5 limiting 
magnitude (Rcfrcgicr et al. 2006), which should detect mas- 
sive spheroids at z ~ 5. The James Webb Space Telescope 
easily has enough sensitivity to detect massive quasar pro- 
genitors out to z > 15 (with stellar masses M* < 1O 1O M0), 
but currently suggested observations, such as a Deep- Wide 
Survey with a 100 arcmin 2 area and limiting m agnitude of 
mAB ~ 30.2 at A = 3.5/xm l|Gardner et al.| [2006). would be 
too narrow to find such rare objects at the highest redshifts. 

6. DISCUSSION 

The discovery of very massive high-redshift galaxies would 
provide an interesting new perspective on the structure for- 
mation process. While the larger abundance of less-massive 
galaxies allow the redshift evolution of such systems to be 
researched by characterizing the statistical properties of pop- 
ulations of objects, the rarity of the most massive galaxies 
will likely limit our ability to draw inferences about their for- 
mation from population studies. If the most massive galaxies 
at the current epoch were a largely heterogeneous population, 
their rarity could prove a significant limitation in unravelling 
the various modes for their formation. 

12 http:/ /pan-starrs. ifa.hawaii.edu/ 

13 http:/ /www. lsst.org/Science/docs/SRD _summary.pdf 

14 http://snap.lbl.gov 



Fortunately, the most massive galaxies today are a roughly 
homogeneous popula tion. Studies of the color bimodal- 
ity of galaxies (e .g., IStrateva etHI 120011 ; iBell et all 120041 ; 
iBaldrv et al.|[2004l ) show that the most luminous galaxies are 
almost uniformly red (g' — r' > 1.2). Such galaxies pop- 
ulate the center of galaxy clusters and their properties are 
roughly uniform between clusters at th e same epoch (e.g., 
iBower et aPl 19921 ; IcTadders fc Yedl2000l ). Importantly, these 
massive galaxies have been demonstrated to contain super- 
massive black holes with mass Mb h > 10 9 M Q (i.e. M87, 
lHarms et aflll994l ; iMacchetto et al.lll997T ). The large SMBH 
masses directly connect the population of the most massive 
galaxies at low redshift to the luminous quasars at z ~ 6, 
since 1O 9 M SMBHs are the most credible engine to power 
2 ~ 6 quasars, and suggest that the most luminous high- 
redshift quasars evolve into a roughly homogeneous popu- 
lation at the current epoch. Our calculations support this 
possibility, under the condition that the most massive galax- 
ies undergo little star formation and their stellar populations 
evolve passively from z ~ 4 to the present. Observed clus- 
ter el lipticals have stellar formation epochs of at least z > 2 
(e.g.. Ivan Dokkum fc Franxll2001l : iGebhardt et al.ll2003l ). We 
note that this picture does not preclude fu ture major mergers 
between sphero ids at lower redshifts fe.g.. Ivan Dokkum 2005; 
IBell et al.ll2006h . which may h ave some relation to core-cusp 
bimodality fe.g.. iFaber et all Il997i : lLauer et all l2006l . Kor- 
mendy et al., in preparation, K rause et al., in preparation, 
though see Fer rarese et al] (|2006T that is the primary hetero- 
geneity in massive galaxies. Even so, this bimodality exists 
primarily in systems with My > —22, with more luminous 
systems mostly displaying surface brightness cores. 

The connection between the most massive galaxies in the 
present epoch and quasars at z ~ 6, forged by their shared 
SMBH masses, and the uniformity of the most massive lo- 
cal galaxies permits an argument for a single evolutionary 
scenario for the formation of the most massive galaxies from 
very high-redshifts to the present day. The appeal of a sin- 
gle mode of formation for the most massive galaxies stems 
primarily from the difficulty in forming luminous quasars at 
high-re dshifts with > 10 9 Mp SM BHs, the problem popular- 
ized by Efstathiou fc Reesl l|1988f ). Finding a robust way of 
growing SMBHs quickly in the limited time available before 
z ~ 6 is problematic, but the success of the ILi et al] (|2006T ) 
simulation of 2 ~ 6 quasar in producing a SMBH with a 
mass of Mbh > 10 9 Mq has provided some hope that high- 
redshift quasars can be explained naturally in the context 
of the formation of rare density peaks in our ACDM cos- 
mology. The calculations performed in this paper provide 
a detailed characterization of the observable ramifications of 
this scenario, the foremost being the possible detection of the 
star-bursting progenitors of 2 ~ 6 quasars at higher redshifts 
(2 > 6) with massive stellar populations (M* ~ irj 11 - 5 - 120 ) 
in wide-area, Lyman-break dropout samples or through wide- 
area IR searches for systems with rest-frame optical breaks at 
2 ~ 9 and stellar masses of Mi, ~ 10 n M Q . 

7. CONCLUSIONS 

Combining hydrodynami cal simulations of the hierarchical 
formation of a 2 ~ 6 quasar (ILi et al . 2006), s tellar population 
synth esis models (ILeitherer et al] (l999; Vaz quez fc Leithererl 
20051). AGN spectr al templates (jVanden Berk et al] 120011 ; 
Marconi et al] [2004) , models for the wavelength-dependent 
attenuation owing to interstellar and inte rgalactic absorption 
l|Calzetti et al.lH99l |2000| ; lMadaulH995T ). and the transmis- 
sivity of telescopes, filters, and detectors, the photometric 
properties of a massive 2 ~ 6 quasar host galaxy are calcu- 
lated at redshifts 2 ~ 4 — 14. The photometric properties of 
the quasar host galaxy reflect three main evolutionary phases 
in its formation; a starburst-dominated phase at z > 6.7, an 
AGN phase at 5.7 < 2 < 6.7 that includes optical quasar 
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activity at z ~ 6.5, and a passive-evolution phase at 2 < 5.7. 
The photometric properties of the system are then compared 
with color-selection techniques for finding high-redshift galax- 
ies. The main findings of these calculations follow. 

• At very high-redshifts z ~ 6.7 — 14, multiple gas-rich 
mergers gives rise to star formation rates of SFR ~ 
10 3 - 10 4 , building a stellar spheroid of > 10 12 M Q by 
z ~ 7 l|Li et al.| [2006) . Our calculations show that these 
star formation rates naturally give rise to a photometric 
starburst phase with strong rest-frame UV and blue op- 
tical emission that satisfy HST ACS/NI CMOS Lyman- 
break criteria for galaxies at z > 7 (e.g. iBouwens et al.l 
120041 ; iBouwens fc Illingworthll2006T ). Over this epoch, 
the near-infrared luminosity of the system peaks at 
z ~ 8.5 with Jixo ~ 21.7, H wo ~ 20.5, and K s ~ 20.5 
in AB magnitudes. 

• During the phase of AGN activity (5.7 < z < 6.7), 
IGM absorption blue-ward of Lya causes the system 
to satisfy both F 775W-band dropout criteria for LBGs 
in GO ODS (e.g. iStanwav et all 120031 ; iGiavalisco et al.1 
l2004aT ) and the SDSS ii'- band dropout s election for 
high-redshift quasars (e.g. ^neTal] [2003). At higher 
redshifts IGM absorption begins to attenuate the z'- 
band flux and push the quasar host galaxy towards the 



locus of L-dwarfs in the z — Jv cg a vs. 



color-color 



space. During this phase, the quasar host gal axy is usu- 
ally ~ 2 z'-band magnitudes fainter than the lFan et al.l 
(2003) sample selection flux limit. 

■ At lower redshifts (2 < 5.7) after the starburst and 
AGN activity has declined, the quasar host galaxy sat- 
isfies V606-dropout crit e ria for galaxies at z ~ 5 (e.g., 
IGiavalisco e t al. 2004tJ lYan et all 12005). The quasar 
descendant also satisfies Balmer break selection tech- 
niques designed t o identify massive and evolved galax- 
ies at 2 ~ 5 (e.g., iWiklind et al1[2007il rbl). 

■ Throughout the photometric evolution of the quasar 
host galaxy, the system satisfies Spitzer IR AC color- 
color selections fo r AGN JStern et al. 2005b) and high- 
redshift galaxies l|Barmbv et a l. 2004). 



The area and depth necessary to find the most mas- 
sive high-redshift galaxies in photometric surveys is 



estimated from our photometric modeling and the 
redshift-dependent abundance of dark matter halos 
(|Press fc Schechter|[l97l ; ISheth fc Torineni ri999). Ex- 
isting surveys, such as the IRAC Shallow Survey or 
NOAO Deep- Wide Field Survey, may be large enough 
to find the massive 2 ~ 4 descendants of 2 ~ 6 quasars. 
Surveys of significant fractions of the sky (> 0.001) in 
the 2850 and Jno or similar bands, such as the UKIRT 
Deep Sky Survey, will be necessary to find the star- 
bursting progenitors of 2 ~ 6 at high-redshifts (2 > 7). 
Future surveys, such as Pan-STARRS, the Large Syn- 
optic Survey Telescope, or the Supernovae Acceleration 
Probe, should have the necessary area and depth to find 
the most massive high-redshift galaxies. 

• If future surveys contain a Jno-dropout massive quasar 
progenitor at 2 ~ 10, evidence for star formation at 
2 > 12 in the form of a Balmer/4000A break could be 
detected in the [3.6pm] — [4.5pm] vs. [4.5pm] — [5.8pm] 
color-color space. 

The remarkable luminosity and early-formation time of the 
massive progenitor galaxies of 2 ~ 6 quasars could enable 
their detection at very high-redshifts 2 > 7. The detection 
of these massive galaxies would help settle outstanding ques- 
tions about how the rarest structures form in the ACDM cos- 
mology and what evolutionary processes give rise to 2 ~ 6 
quasars. As observations must ultimately determine the exis- 
tence and properties of the most massive high-redshift galax- 
ies, we eagerly anticipate future efforts to detect and charac- 
terize these extraordinary quasar progenitors. 



BER gratefully acknowledges support from a Spitzer Fel- 
lowship, through a NASA grant administered by the Spitzer 
Science Center, and thanks Hsiao- Wen Chen, Tiziana Di Mat- 
teo, Volker Springel, Marijn Franx, and Andrew Zentner for 
helpful comments. BER also thanks Dr. Wiklind for early 
access to his results before publication. This work was sup- 
ported in part by NSF grants ACI 96-19019, AST 00-71019, 
AST 02-06299, and AST 03-07690, and NASA ATP grants 
NAG5-12140, NAG5-13292, and NAG5-13381. The simula- 
tions were performed at the Center for Parallel Astrophysical 
Computing at Harvard-Smithsonian Center for Astrophysics. 



REFERENCES 



Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 
Baldry, I. K., Glazebrook, K., Brinkmann, J., Ivezic, Z., Lupton, 

R. H., Nichol, R. C, & Szalay, A. S. 2004, ApJ, 600, 681 
Barmby, P., et al. 2004, ApJS, 154, 97 
Barnes, J. E., & Hernquist, L. 1996, ApJ, 471, 115 
Barnes, J. E., & Hernquist, L. E. 1991, ApJ, 370, L65 
Barth, A. J., Ho, L. C, & Sargent, W. L. W. 2003, ApJ, 583, 134 
Barton, E. J., Dave, R., Smith, J.-D. T., Papovich, C, Hernquist, 

L., & Springel, V. 2004, ApJ, 604, LI 
Baugh, C. M., Cole, S., Frenk, C. S., & Lacey, C. G. 1998, ApJ, 

498, 504 

Beckwith, S. V. W., et al. 2006, AJ, 132, 1729 
Bell, E. P., et al. 2006, ApJ, 640, 241 
— . 2004, ApJ, 608, 752 

Bershady, M. A., Charlton, J. C, & Geoffroy, J. M. 1999, ApJ, 
518, 103 

Bertoldi, F., Carilli, C. L., Cox, P., Fan, X., Strauss, M. A., Beelen, 

A., Omont, A., & Zylka, R. 2003, A&A, 406, L55 
Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 
Bondi, H. 1952, MNRAS, 112, 195 

Bouwens, R. J., & Illingworth, G. D. 2006, Nature, 443, 189 
Bouwens, R. J., et al. 2003, ApJ, 595, 589 
— . 2004, ApJ, 616, L79 

Bower, R. G., Lucey, J. R., & Ellis, R. S. 1992, MNRAS, 254, 601 
Bremer, M. N., Jensen, J. B., Lehnert, M. D., Schreiber, N. M. F., 

& Douglas, L. 2004, ApJ, 615, LI 
Brodwin, M., et al. 2006, ApJ, 651, 791 



Bromm, V., k. Larson, R. B. 2004, ARA&A, 42, 79 

Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., Kravtsov, 

A. V., Klypin, A. A., Primack, J. R., & Dekel, A. 2001, MNRAS, 
321, 559 

Bunker, A. J., Stanway, E. R., Ellis, R. S., & McMahon, R. G. 

2004, MNRAS, 355, 374 
Bunker, A. J., Stanway, E. R., Ellis, R. S., McMahon, R. G., k, 

McCarthy, P. J. 2003, MNRAS, 342, L47 
Calzetti, D. 1997, AJ, 113, 162 

Calzctti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornneef, 
J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 
429, 582 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Carilli, C. L., et al. 2004, AJ, 128, 997 

Chary, R.-R., Stern, D., & Eisenhardt, P. 2005, ApJ, 635, L5 
Chen, H.-W., & Marzke, R. O. 2004, ApJ, 615, 603 
Cohen, M., Wheaton, W. A., & Megeath, S. T. 2003, AJ, 126, 1090 
Coil, A. L., Newman, J. A., Cooper, M. C, Davis, M., Faber, S. M., 

Koo, D. C, & Willmer, C. N. A. 2006, ApJ, 644, 671 
Cowie, L. L., & Hu, E. M. 1998, AJ, 115, 1319 
Cox, A. N. 2000, Allen's astrophysical quantities (Allen's 

astrophysical quantities, 4th ed. Publisher: New York: AIP Press; 

Springer, 2000. Editedy by Arthur N. Cox. ISBN: 0387987460) 
Cox, T. J., Di Matteo, T., Hernquist, L., Hopkins, P. F., Robertson, 

B. , & Springel, V. 2006a, ApJ, 643, 692 



18 



Robertson et al. 



Cox, T. J., Dutta, S. N., Di Matteo, T., Hernquist, L., Hopkins, 

P. F., Robertson, B., & Springel, V. 2006b, ApJ, 650, 791 
Croom, S. M., et al. 2005, MNRAS, 356, 415 
Daddi, E., et al. 2003, ApJ, 588, 50 

Dave, R., Finlator, K., & Oppenheimer, B. D. 2006, MNRAS, 370, 
273 

Dave, R., Hernquist, L., Katz, N., & Weinberg, D. H. 1999, ApJ, 
511, 521 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

Dickinson, M., Giavalisco, M., & The Goods Team. 2003, in The 
Mass of Galaxies at Low and High Redshift, ed. R. Bender & 

A. Renzini, 324-+ 

Dickinson, M., et al. 2004, ApJ, 600, L99 

Dunlop, J. S., Cirasuolo, M., & McLurc, R. J. 2006, ArXiv 

Astrophysics e-prints 
Efstathiou, G., & Rees, M. J. 1988, MNRAS, 230, 5P 
Egami, E., et al. 2005, ApJ, 618, L5 
Eisenhardt, P. R., et al. 2004, ApJS, 154, 48 

Ellis, R., Santos, M. R., Kneib, J.-P., & Kuijken, K. 2001, ApJ, 
560, L119 

Eyles, L., Bunker, A., Ellis, R., Lacy, M., Stanway, E., Stark, D., 

& Chiu, K. 2006, ArXiv Astrophysics c-prints 
Eyles, L. P., Bunker, A. J., Stanway, E. R., Lacy, M., Ellis, R. S., 

& Doherty, M. 2005, MNRAS, 364, 443 
Faber, S. M., et al. 1997, AJ, 114, 1771 
Fabian, A. C. 1999, MNRAS, 308, L39 

Fagotto, F., Bressan, A., Bertelli, G., & Chiosi, C. 1994, A&AS, 
105, 29 

Fan, X., et al. 2004, AJ, 128, 515 

— . 2001a, AJ, 122, 2833 

— . 2006, AJ, 131, 1203 

— . 2001b, AJ, 121, 31 

— . 2003, AJ, 125, 1649 

— . 2000, AJ, 120, 1167 

Fauchcr-Gigucre, C.-A., Lidz, A., Zaldarraiga, M., & Hernquist, L. 

2007, ArXiv Astrophysics e-prints 
Fazio, G. G., ct al. 2004, ApJS, 154, 10 
Ferrarese, L., et al. 2006, ApJS, 164, 334 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Finlator, K., Dave, R., Papovich, C, & Hernquist, L. 2006, ApJ, 

639, 672 

Fitzpatrick, E. L. 1986, AJ, 92, 1068 

Ford, H. C, et al. 1998, in Proc. SPIE Vol. 3356, p. 234-248, 
Space Telescopes and Instruments V, ed. P. Y. Bely & J. B. 
Breckinridge 

Franx, M., et al. 2003, ApJ, 587, L79 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., 

& Schneider, D. P. 1996, AJ, 111, 1748 
Gao, L., White, S. D. M., Jenkins, A., Frenk, C. S., & Springel, V. 

2005, MNRAS, 363, 379 
Gardner, J. P., et al. 2006, Space Science Reviews, 123, 485 
Gebhardt, K., ct al. 2000, ApJ, 539, L13 
— . 2003, ApJ, 597, 239 
Genzel, R., et al. 1998, ApJ, 498, 579 
Giavalisco, M., et al. 2004a, ApJ, 600, L103 
— . 2004b, ApJ, 600, L93 

Gladders, M. D., & Yee, H. K. C. 2000, AJ, 120, 2148 

Gordon, K. D., Clayton, G. C, Misselt, K. A., Landolt, A. U., & 

Wolff, M. J. 2003, ApJ, 594, 279 
Granato, G. L., De Zotti, G., Silva, L., Bressan, A., & Danese, L. 

2004, ApJ, 600, 580 
Grazian, A., et al. 2006, A&A, 449, 951 

Guimaraes, R. N., Petitjean, P., Rollinc, E. B., De Carvalho, R. R., 
Djorgovski, G., Srianand, R., Aghaee, A., & Castro, S. 2007, 
ArXiv Astrophysics e-prints 

Harford, A. G., & Gnedin, N. Y. 2003, ApJ, 597, 74 

Haring, N., & Rix, H.-W. 2004, ApJ, 604, L89 

Harms, R. J., et al. 1994, ApJ, 435, L35 

Hernquist, L. 1989, Nature, 340, 687 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B. , & Springel, V. 2006a, ApJS, 163, 1 

Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Krausc, 

E. 2007a, ArXiv Astrophysics e-prints 
Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & Springel, 

V. 2006b, ApJS, 163, 50 
Hopkins, P. F., Hernquist, L., Martini, P., Cox, T. J., Robertson, 

B., Di Matteo, T., & Springel, V. 2005, ApJ, 625, L71 



Hopkins, P. F., Lidz, A., Hernquist, L., Coil, A. L., Myers, A. D., 
Cox, T. J., & Spergel, D. N. 2006c, ArXiv Astrophysics c-prints 

Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007b, ApJ, 654, 
731 

Hu, E. M., & Cowie, L. L. 2006, Nature, 440, 1145 

Hu, E. M., Cowie, L. L., McMahon, R. G., Capak, P., Iwamuro, 

F., Kneib, J.-P., Maihara, T., & Motohara, K. 2002, ApJ, 568, 

L75 

Irwin, M., McMahon, R. G., & Hazard, C. 1991, in ASP Conf. 
Ser. 21: The Space Distribution of Quasars, ed. D. Crampton, 
117-126 

lye, M., et al. 2006, Nature, 443, 186 

Jannuzi, B. T., & Dey, A. 1999, in ASP Conf. Ser. 191: Photometric 
Redshifts and the Detection of High Redshift Galaxies, ed. 
R. Weymann, L. Storrie-Lombardi, M. Sawicki, & R. Brunner, 
111-+ 

Kashikawa, N., et al. 2006, ApJ, 648, 7 

Katz, N., Hernquist, L., & Weinberg, D. H. 1992, ApJ, 399, L109 
Kauffmann, G., Colberg, J. M., Diaferio, A., & White, S. D. M. 

1999, MNRAS, 303, 188 
Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 311, 576 
Kennicutt, R. C. 1998, ApJ, 498, 541 

Kim, Y.-R., & Croft, R. 2007, ArXiv Astrophysics e-prints 

King, A. 2003, ApJ, 596, L27 

Kodaira, K., et al. 2003, PASJ, 55, L17 

Koester, B. P., et al. 2007, ArXiv Astrophysics e-prints 

Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 

Kroupa, P. 2001, MNRAS, 322, 231 

— . 2002, Science, 295, 82 

Kurk, J. D., Cimatti, A., di Serego Alighieri, S., Vernet, J., Daddi, 
E., Fcrrara, A., & Ciardi, B. 2004, A&A, 422, L13 

Labbe, I., Bouwens, R., Illingworth, G. D., & Franx, M. 2006, ApJ, 
649, L67 

Labbe, I., et al. 2005, ApJ, 624, L81 

Lanzetta, K. M., Yahil, A., & Fernandez-Soto, A. 1996, Nature, 
381, 759 

Larson, R. B., & Tinsley, B. M. 1978, ApJ, 219, 46 

Lauer, T. R., et al. 2006, ArXiv Astrophysics e-prints 

Lawrence, A., et al. 2006, ArXiv Astrophysics e-prints 

Lee, K.-S., Giavalisco, M., Gnedin, O. Y., Somerville, R. S., 

Ferguson, H. C, Dickinson, M., & Ouchi, M. 2006, ApJ, 642, 

63 

Lehnert, M. D., & Bremer, M. 2003, ApJ, 593, 630 
Leitherer, C, et al. 1999, ApJS, 123, 3 
Li, Y., ct al. 2006, ArXiv Astrophysics e-prints 
Lynden-Bell, D. 1969, Nature, 223, 690 

Macchetto, F., Marconi, A., Axon, D. J., Capctti, A., Sparks, W., 

& Crane, P. 1997, ApJ, 489, 579 
Madau, P. 1995, ApJ, 441, 18 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., 

Steidel, C. C, & Fruchter, A. 1996, MNRAS, 283, 1388 
Magdziarz, P., & Zdziarski, A. A. 1995, MNRAS, 273, 837 
Magorrian, J., et al. 1998, AJ, 115, 2285 

Maiolino, R., Schneider, R., Oliva, E., Bianchi, S., Fcrrara, A., 
Mannucci, F., Pedani, M., & Roca Sogorb, M. 2004, Nature, 
431, 533 

Marconi, A., et al. 2003, ApJ, 586, 868 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., & 

Salvati, M. 2004, MNRAS, 351, 169 
McKee, C. F., & Ostriker, J. P. 1977, ApJ, 218, 148 
Mihos, J. C, & Hernquist, L. 1994, ApJ, 431, L9 
— . 1996, ApJ, 464, 641 

Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Mobasher, B., et al. 2005, ApJ, 635, 832 

Nagaminc, K., Ccn, R., Furlanetto, S. R., Hernquist, L., Night, C, 
Ostriker, J. P., & Ouchi, M. 2006, New Astronomy Review, 50, 
29 

Nagaminc, K., Cen, R., Hernquist, L., Ostriker, J. P., & Springel, 

V. 2005a, ApJ, 627, 608 
— . 2005b, ApJ, 618, 23 

Nagaminc, K., Springel, V., & Hernquist, L. 2004, MNRAS, 348, 
421 

Narayanan, D., ct al. 2006, ApJ, Submitted 

Night, C, Nagamine, K., Springel, V., & Hernquist, L. 2006, 

MNRAS, 366, 705 
Noguchi, M. 1988, A&A, 203, 259 
Oke, J. B. 1974, ApJS, 27, 21 
Ouchi, M., et al. 2005, ApJ, 620, LI 



Massive High-Redshift Galaxies 



19 



— . 2004, ApJ, 611, 660 

Overzier, R. A., Bouwens, R. J., Illingworth, G. D., & Pranx, M. 

2006, ArXiv Astrophysics e-prints 
Papovich, C, Dickinson, M., & Ferguson, H. C. 2001, ApJ, 559, 

620 

Papovich, C, et al. 2006, ApJ, 640, 92 
Pci, Y. C. 1992, ApJ, 395, 130 

Porciani, C, Magliocchetti, M., & Norberg, P. 2004, MNRAS, 355, 
1010 

Porciani, C, & Norberg, P. 2006, MNRAS, 371, 1824 

Power, C, Navarro, J. F., Jenkins, A., Frcnk, C. S., White, 

S. D. M., Springel, V., Stadel, J., & Quinn, T. 2003, MNRAS, 

338, 14 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Reddy, N. A., & Steidel, C. C. 2004, ApJ, 603, L13 

Reddy, N. A., Steidel, C. C, Erb, D. K., Shaplcy, A. E., & Pettini, 

M. 2006, ArXiv Astrophysics c-prints 
Refregier, A., et al. 2006, in Proc. SPIE Vol. 6265, p. 62651Y, 

Space Telescopes and Instrumentation I, ed. J. C. Mather, H. A. 

MacEwen & M. W. M. de Graauw 
Richard, J., Pello, R., Schaerer, D., Lc Borgnc, J.-F., & Kneib, 

J. -P. 2006, ArXiv Astrophysics e-prints 
Robertson, B., Cox, T. J., Hernquist, L., Franx, M., Hopkins, P. F., 

Martini, P., & Springel, V. 2006a, ApJ, 641, 21 
Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, T., Hopkins, 

P. F., Martini, P., & Springel, V. 2006b, ApJ, 641, 90 
Rodrigues, C. V., Magalhaes, A. M., Coyne, G. V., & Piirola, V. 

1997, ApJ, 485, 618 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 
Sanders, D. B., Soifer, B. T., Elias, J. H., Madore, B. F., Matthews, 

K., Neugebauer, G., & Scoville, N. Z. 1988a, ApJ, 325, 74 
Sanders, D. B., Soifer, B. T., Elias, J. H., Neugebauer, G., & 

Matthews, K. 1988b, ApJ, 328, L35 
Santos, M. R. 2004, MNRAS, 349, 1137 
Schmidt, M. 1959, ApJ, 129, 243 
Schneider, D. P., et al. 2001, AJ, 121, 1232 

Schneider, D. P., Schmidt, M., & Gunn, J. E. 1991, AJ, 101, 2004 
Shapley, A. E., Steidel, C. C, Erb, D. K., Reddy, N. A., 

Adelberger, K. L., Pettini, M., Barmby, P., & Huang, J. 2005, 

ApJ, 626, 698 
Shcn, Y., et al. 2007, ArXiv Astrophysics e-prints 
Shcth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 
Silk, J., & Rees, M. J. 1998, A&A, 331, LI 

Skrutskie, M. F., ct al. 1997, in ASSL Vol. 210: The Impact of 
Large Scale Near-IR Sky Surveys, ed. F. Garzon, N. Epchtcin, 
A. Omont, B. Burton, & P. Persi, 25-+ 

Soifer, B. T., et al. 1984, ApJ, 278, L71 

Soifer, B. T., Sanders, D. B., Madore, B. F., Neugebauer, G., 
Daniclson, G. E., Elias, J. H., Lonsdale, C. J., & Rice, W. L. 
1987, ApJ, 320, 238 

Somervillc, R. S., Primack, J. R., & Fabcr, S. M. 2001, MNRAS, 
320, 504 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

Spinrad, H., Stern, D., Bunker, A., Dey, A., Lanzetta, K., Yahil, 

A., Pascarelle, S., & Fernandez-Soto, A. 1998, AJ, 116, 2617 
Springel, V. 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005a, ApJ, 620, L79 
— . 2005b, MNRAS, 361, 776 

Springel, V., & Hernquist, L. 2003, MNRAS, 339, 289 
Stanway, E. R., Bunker, A. J., & McMahon, R. G. 2003, MNRAS, 
342, 439 

Stanway, E. R., Bunker, A. J., McMahon, R. G., Ellis, R. S., Treu, 

T., & McCarthy, P. J. 2004a, ApJ, 607, 704 
Stanway, E. R., ct al. 2004b, ApJ, 604, L13 

Stanway, E. R., McMahon, R. G., & Bunker, A. J. 2005, MNRAS, 
359, 1184 

Stark, D. P., Bunker, A. J., Ellis, R. S., Eylcs, L. P., & Lacy, M. 
2006, ArXiv Astrophysics e-prints 



Steidel, C. C, Pettini, M., & Hamilton, D. 1995, AJ, 110, 2519 
Steidel, C. C, Giavalisco, M., Pettini, M., Dickinson, M., & 

Adelberger, K. L. 1996, ApJ, 462, L17+ 
Steidel, C. C, & Hamilton, D. 1992, AJ, 104, 941 
— . 1993, AJ, 105, 2017 
Stern, D., et al. 2005a, ApJ, 631, 163 

Stern, D., Yost, S. A., Eckart, M. E., Harrison, F. A., Helfand, 
D. J., Djorgovski, S. G., Malhotra, S., & Rhoads, J. E. 2005b, 
ApJ, 619, 12 

Stiavelli, M., Fall, S. M., & Panagia, N. 2004, ApJ, 610, LI 
Strateva, I., et al. 2001, AJ, 122, 1861 
Taniguchi, Y., et al. 2005, PASJ, 57, 165 

Telfer, R. C, Zheng, W., Kriss, G. A., & Davidsen, A. F. 2002, 

ApJ, 565, 773 
Thompson, R. I., et al. 2005, AJ, 130, 1 

Toomrc, A. 1977, in Evolution of Galaxies and Stellar Populations, 
401— I- 

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 

Tremaine, S., et al. 2002, ApJ, 574, 740 

van Dokkum, P. G. 2005, AJ, 130, 2647 

van Dokkum, P. G., & Franx, M. 2001, ApJ, 553, 90 

van Dokkum, P. G., et al. 2004, ApJ, 611, 703 

Vanden Berk, D. E., et al. 2001, AJ, 122, 549 

Vanzella, E., et al. 2006, A&A, 454, 423 

Vazquez, G. A., & Leitherer, C. 2005, ApJ, 621, 695 

Verma, A., Lehnert, M. D., Foerster Schreiber, N. M., Bremer, 

M. N., & Douglas, L. 2007, ArXiv Astrophysics e-prints 
Vignali, C, Brandt, W. N., & Schneider, D. P. 2003, AJ, 125, 433 
Volonteri, M., Haardt, F., & Madau, P. 2003, ApJ, 582, 559 
Wake, D. A., et al. 2004, ApJ, 610, L85 
Walter, F., et al. 2003, Nature, 424, 406 

Walter, F., Carilli, C, Bertoldi, F., Menten, K., Cox, P., Lo, K. Y., 

Fan, X., & Strauss, M. A. 2004, ApJ, 615, L17 
Warren, S., & Hewett, P. 2002, in ASP Conf. Ser. 283: A New Era 

in Cosmology, ed. N. Metcalfe & T. Shanks, 369-+ 
Warren, S. J., et al. 2007, MNRAS, 375, 213 

Warren, S. J., Hewett, P. C, Irwin, M. J., McMahon, R. G., & 

Bridgeland, M. T. 1987, Nature, 325, 131 
Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, A. V., & 

Dekel, A. 2002, ApJ, 568, 52 
Weinberg, D. H., Hernquist, L., & Katz, N. 2002, ApJ, 571, 15 
Werner, M. W., et al. 2004, ApJS, 154, 1 

Wiklind, T., Mobasher, B., Dickinson, M., Ferguson, H. C, 

Giavalisco, M., Grogin, N. A., & Panagia, N. 2007a, IAU 

Symposium, 235, 368 
Wiklind, T., Dickinson, M., Ferguson, H. C, Giavalisco, M., 

Mobasher, B., Grogin, N. A., & Panagia, N. 2007b, ApJ, 

submitted. 

Williams, R. E., et al. 1996, AJ, 112, 1335 

Willott, C. J., Rawlings, S., Jarvis, M. J., & Blundell, K. M. 2003, 

MNRAS, 339, 173 
Wyithe, J. S. B., & Loeb, A. 2003, ApJ, 595, 614 
Yan, H., et al. 2004, ApJ, 616, 63 

Yan, H., Dickinson, M., Giavalisco, M., Stern, D., Eisenhardt, 

P. R. M., & Ferguson, H. C. 2006, ApJ, 651, 24 
Yan, H., ct al. 2005, ApJ, 634, 109 
Yan, H., & Windhorst, R. A. 2004, ApJ, 612, L93 
Yan, H., Windhorst, R. A., & Cohen, S. H. 2003, ApJ, 585, L93 
York, D. G., ct al. 2000, AJ, 120, 1579 

Yoshida, N., Abel, T., Hernquist, L., & Sugiyama, N. 2003, ApJ, 
592, 645 

Yoshida, N., Omukai, K., Hernquist, L., & Abel, T. 2006, ApJ, 652, 
6 

Young, J. S., Schlocrb, F. P., Kenney, J. D., & Lord, S. D. 1986, 

ApJ, 304, 443 
Zamorani, G., et al. 1981, ApJ, 245, 357 



